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Research  Objective:  The  objective  of  this  research  program  is  to  develop  a  fundamental 
understanding  of  the  gas-phase  reaction  mechanisms  associated  with  the  combustion  pro¬ 
cess  of  energetic  polycyclic  mtramines  and  energetic  polymers. 

Summar>’  of  .Major  Accomplishments: 


1.  Kinetic  experiments  on  N^O  and  NO2  decomposition  and  on  the  reactions  betw'een 
N2O  and  NO2  with  H2.  all  with  dilute  mixmres  in  N2.  have  been  performed  in  a  flow 
reactor  over  the  temperature  range  of  803-1100  K  and  pressure  range  of  1-6  atm.  Spe¬ 
cies  concentration  profiles  measured  include  H2.  H2O.  NO.  NO2.  N2O.  and  O2. 

2.  A  detailed  chemical  kinetic  model  has  been  constructed  from  literature  rate  constants 
and  model  predictions  have  been  compared  with  experimental  data.  Agreement  be¬ 
tween  model  predictions  and  experiment  for  the  H2/N2O  system  were  generally  good, 
while  for  the  H2/NO2  system,  the  model  predicts  a  reaction  rate  approximately  a  factor 
of  2.5  times  too  fast. 

3.  Sensitivity  and  reaction  flux  analyses  have  been  performed  and  crucial  reactions  requir¬ 
ing  further  theoretical  and  experimental  study  have  been  identified. 

4.  For  the  H2/N2O  system,  the  reaction  OH  -1-  N2O  ->  HO2  +  N2  was  found  to  be  an  im¬ 
portant  reaction  leading  to  inhibition.  This  reaction  produces  HO2  radicals  which 
makes  their  kinetics  also  important  to  the  system.  The  channel  H  +  N2O  ->  NH  +  NO 
was  also  observed  to  inhibit  the  reaction,  but  to  have  a  much  smaller  effect.  Nitric  oxide 
formation  was  not  observed  in  the  experimental  data  and  hence  this  latter  step  may  not 
be  important  for  the  present  conditions. 


5.  For  the  H2/NO2  system.  HONO  formation  and  consumption  were  found  to  be  impor¬ 
tant.  The  rates  of  the  initiation  reaction  H2  +  NO2  ->  HONO  +  H  and  the  reaction  • 

HONO  +  OH  ->  H2O  +  NO2  are  likely  responsible  for  the  observed  discrepancies  be¬ 
tween  model  and  experiment. 


Abstract 


The  objective  of  this  research  program  is  to  develop  a  fundamental  understanding  of  the 
gas-phase  reaction  mechanisms  associated  with  the  combustion  process  of  energetic  poly¬ 
cyclic  nitramines  and  energetic  polymers.  Kinetic  experiments  are  performed  in  a  variable 
temperature  (600  K  <T<  1200  K),  variable  pressure  (1  ami.  <  P  <  18  ami.)  flow  reactor 
with  residence  times  ranging  from  10  mi  to  10  s.  Reaction  mechanisms,  based  on  the  prin¬ 
ciple  of  hierarchical  model  construction,  are  developed  and  validated  by  detailed  kinetic 
modeling  of  the  expenmenis  along  with  reaction  flux  and  sensitivity  analyses.  During  the 
past  year,  a  Nicolet  Model  730  FTiR  was  installed  and  calibrated  for  the  measurement  of 
oxides  of  nitrogen,  simple  hydrocarbons.  CO,  CO2,  and  H2O.  Under  other  suppon.  the  de¬ 
velopment  of  a  continuous  analytical  techruque  for  hydrogen  was  also  developed.  Experi¬ 
ments  and  modeling  calculations  were  conducted  on  the  kinetics  of  N2O  and  NO2 
decomposition  (1000  K  -1100  K.  3  ami.  -  6  ami.),  and  on  the  reactions  of  H2  with  N2O 
(925  K  -  1073  K,  1  ami.  -  3  ami.)  and  H2  with  NO2  (803  K  -  873  K.  3  ami.).  These  experi¬ 
ments  and  models  represent  the  basis  for  current  and  future  gas-phase  kinetic  studies  on 
the  chemical  systems  consisting  of  the  fuels  CO,  CH2O.  C//4.  C2//4.  HCN.  and  NH-^  and 
the  oxidizers  NO2  and  N2O. 
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1.0  INTRODUCTION 


1.1  Statement  of  the  Problem  and  Research  Objectives 

Newly  developed  energetic  polycyclic  nitramines  and  energetic  polymers  have  the  poten¬ 
tial  to  offer  significant  energy  density  and  plume  signature  advantages  over  current  state- 
of-the-an  propellants  in  strategic,  tactical  rocket  and  gun  propulsion.  Current  state-of-the- 
art  propellants,  such  as  HMX  and  RDX,  were  developed  empirically  over  a  lengthy  period 
of  time,  being  first  synthesized  nearly  a  century  ago.  The  development  and  implementa¬ 
tion  of  the  next  generation  of  propellants  will  have  to  meet  new  and  more  stringent  re¬ 
quirements  of  energy  density,  bum  rate,  stability,  signature  and  hazards  specifications,  and 
formulations  will  likely  be  based  on  combinations  of  more  than  ten  new  oxidizers  and  en¬ 
ergetic  polymers.  To  provide  some  focal  direction  in  what  will  be  a  lengthy  (and  otherwise 
nsky)  empirical  development,  a  national  energetic  materials  research  program  is  emerg¬ 
ing. 

The  goal  of  this  national  research  program  is  to  develop  a  robust,  comprehensive  chemical 
and  physical  model  for  the  combustion  of  new  polycyclic  nitramines.  The  proposed  pro¬ 
gram  will  provide  a  tool  which  can  be  used  not  only  for  interpreting  experimental  results, 
but  for  predicting  quantitatively  the  effects  of  changing  propellant  design.  Potential  bene¬ 
fits  include  reductions,  both  in  expensive  laboratory  and  full  scale  testing,  and  in  the  time 
required  for  successful  development.  The  development  of  such  a  model  will  require  a  fun¬ 
damental  quantitative  understanding  of  the  elementary  chemical  and  physical  processes  in 
the  gas.  liquid  and  solid  phases  of  propellant  ignition  and  combustion.  The  currently  pro¬ 
posed  approach  for  development  of  this  comprehensive  model  is  based  upon  a  hierarchical 
structure  as  shown  in  FIGURE  1  and  consists  of  submodel  testing  and  validation  and  the  ap¬ 
plication  of  imponant  new  experimental  and  theoretical  tools  developed  for  and  success¬ 
fully  used  in  other  complex  physicochemical  problems.  Among  these  tools  are  state-of- 
the-art  expenmental  techmques  in  chemical  physics  and  physical  chemistry  (e.g..  laser 
spectroscopic  techniques),  state-of-ihe-an  computational  techniques  in  fluid  mechanics, 
chemistry,  and  physics  (largely  due  to  the  development  of  massively  parallel  computer  ar¬ 
chitecture  and  algorithms),  and  systematic  sensitivity  analysis  techniques,  which  link  ex¬ 
perimental  and  computational  microscopic  beha  lor  to  macroscopic  observations. 

It  is  important  to  recognize  that  the  combustion  of  high  energy  density  propellants  is  an¬ 
other  example  of  complex  physicochernical  processes  which  have  demanded  the  research 
community's  attention  in  recent  years.  The  combustion  of  hydrocarbons  and  the  formation 
of  atmospheric  air  pollution  are  two  examples  of  complex  phenomena  where  analog  pro¬ 
grams  have  contributed  immensely  to  the  understanding  and  control  of  'uch  processes. 
The  work  of  Melius  (1988)  on  RDX  combustion  and  of  Yetter  et  at.  (1988)  on  boron  parti¬ 
cle  combustion  are  two  examples  of  detailed  modeling/sensitivity  analysis  studies  on  pro¬ 
pellant  combustion.  With  the  transfer  of  knowledge  from  analog  programs,  the  present 
program  has  a  high  probability  of  success. 

A  fundamental  input  to  the  overall  program  is  the  gas-phase  nitrogen  combustion  chemis¬ 
try  and  transpoit  network  simulation  and  validation  subprogram  at  the  base  of  FlGlTE  1. 
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FIGURE  1.  High  Performance  Propellant  Combustion  Simulation 
Research  Program 


The  goal  of  this  subprogram  is  to  develop  a  detailed  chemical/physical  model  for  reaction 
and  diffusion  of  the  gas-phase  species  associated  with  the  combustion  process  of  energetic 
polycyclic  mcramines  and  energeuc  polymers.  Like  the  overall  program,  this  subprogram 
can  be  subdivided  into  the  hierarchy  as  shown  in  the  flow  chan  of  FIGURE  2.  The  proposed 
approach  stans  with  developing  a  basic  understanding  of  the  gas-  phase  thermochemical, 
kinetic,  and  transpon  properties  of  key  nitramine  species  and  ends  with  the  validation  of  a 
comprehensive  reaction  mechanism  for  the  gas  phase  reaction  of  intermeuiate  species.  In 
this  hierarchical  approach,  theory  at  the  most  fundamental  level  and  at  a  semi-empirical 
level,  interact  with  laboratory  measurements  of  rate  constants,  product  pathways  and 
transport  properties.  From  these  pieces  of  information,  a  well  founded  chemical  mecha¬ 
nism  with  associated  rate  constants  is  first  derived  and  validated  by  comparison  of  model 


Necessary  and  Sufficient 
Experiments 


Coupled  Combustion 
Chemistry,  Transport,  and 

Heat  Rolease  Expehments 

II 

■ 

Diffusion  Flames 

Premixed  Flames 

t 

Combustion  Chemistry 
Expenments 

Shock  Tubes 

Flow  Reactors 

•l 

Isola'ed  Species  and 
Isolaied  Reaction 
Experiments 


Rale  Constant  Measurements 
Product  Detection 
Transport  (Radicals) 


Mjftipie  and  Complimenia^ 
Diagnostics 


Model/Theory 


1-D  Gas  Phase  Model 
Simulation  and  Validation 


Transport  and  Heat  Release 
added  to  Elementary 
Chemistry 


0-D  Gas  Phase  Model 
Simulation  and  Validation 


Complex  Chemistry  Only 


Semi-  Er, pineal  Theory 


Rate  Constant  Estimates 
Cnticat  Data  Analysis 
Thermochemistry 


Electronic  Structure 
Reaction  Rate 
Transport 


Sensitivity  Analys  s 


FIGLRE  2.  Nitrogen  Combustion  Chemistry  and  Transport 
Network  Simulation  and  N'alidation 


prediction  with  kinetic  data  from  shock  tubes  and  flow  reactors.  This  accomplished,  the 
effects  of  heat  release  and  species/energy  transport  are  included  by  simulation  and  valida¬ 
tion  of  premixed  and  diffusion  flames.  Sensitivity  analysis  techniques  are  used  at  each  lev¬ 
el  of  theory  and  model  development  and  validation.  In  panicular,  sensitivity  analysis 
procedures  are  used  to  identify  key  species  and  reactions,  and  hence,  provide  an  interac¬ 
tive  and  guiding  tool  between  the  elemental  and  mechanistic  kinetic  studies.  In  addition, 
multiple  and  complementary  diagnostic  techniques  are  applied  to  each  of  the  experimental 
programs.  The  specific  tasks  of  this  subprogram  include: 

*  determine  via  fundamental  and  semi-empincal  theory  the  thermochemistry  and  trans¬ 
port  properties  of  key  species  and  the  rate  constants  of  key  reactions. 


•  experimentally  measure  rate  constants  and  reaction  products  of  key  reactions  and  trans- 
pon  coefficients  of  key  species. 

•  perform  critical  review  and  analysis  of  therm oc hem ical.  kinetic  and  transport  data  of 
key  species  and  reactions  to  yield  well-founded  reaction  mechamsms. 

•  validate  reaction  mechanisms  and  associated  rate  parameters  over  well-defined  and 
controllable  temperature  and  pressure  ranges  with  transport-free  kinetic  data  from 
shock  tube  experiments  and  flow  reactor  experiments. 

•  validate  combustion  chemistry  with  simultaneous  large  heat  release  with  data  from  gas¬ 
eous  premixed  flames. 

•  validate  transpon  properties  with  data  from  gaseous  counter-flow  diffusion  flames. 

The  assembly  and  validation  of  detailed  kinetic  mechanisms  is  an  important  aspect  of  de¬ 
veloping  understanding  of  chemical  processes  occurring  in  propellant  combustion.  The  re¬ 
sulting  mechamsms  provide  useful  tools  in  determining  what  elementary  reacuons  require 
improved  definition,  in  acting  as  benchmarks  against  which  lumped  chemistry  models  can 
be  developed  and  tested,  and  in  evaluating  the  interactions  of  chemistry  and  transpon  phe¬ 
nomena  in  simple  one-dimensional  systems.  Moreover,  comprehensive  mechamsms.  eval¬ 
uated  and  validated  over  larger  ranges  in  pressure,  temperature,  and  equivalence  ratio  are 
required  for  simulating  the  more  practical  combustion  environments  where  this  gas  phase 
chemistry  is  intimately  coupled  with  chemical  processes  occurring  at  the  material  surface. 

In  the  analog  field  of  hydrocarbon  combustion,  the  oxidation  reaction  consists  primarily  of 
the  sequential  fragmentation  of  the  initial  fuel  molecule  into  smaller  intermediate  species, 
which  are  ultimately  converted  to  final  products,  usually  dominated  by  H2O  and  CO2.  In 
many  cases,  these  intermediate  species  can  be  treated  as  fuels  themselves.  For  example, 
ethylene  is  an  imponant  intermediate  of  propane  combustion  and  other  higher  hydrocar¬ 
bons.  but  ethylene  can  also  be  a  primary  fuel.  Carbon  monoxide  and  hydrogen  are  com¬ 
mon  species  which  are  observed  during  the  oxidation  of  all  hydrocarbons,  and  the  same 
radical  species  H,  O.  OH.  HO2,  HCO.  and  others  are  common  to  all  hydrocarbon  combus¬ 
tion.  These  characteristics  can  be  used  to  great  advantage  in  constructing  reaction  mecha¬ 
nisms  for  complex  practical  fuels.  Reacuon  mechanisms  can  be  developed  systematically, 
beginning  with  the  simplest  species  and  reactions  which  are  common  sub-elements  in  the 
combustion  of  more  complex  species,  and  sequentially  constructed  by  incorporating  new 
species  and  reactions  in  order  of  increasing  complexity.  At  each  level,  the  newly  added 
portions  of  the  mechanism  must  be  tested  and  validated  by  thorough  comparison  between 
numerically  predicted  and  expenmentaJly  observed  results.  However,  because  of  the  se¬ 
quential  ordering,  only  those  features  which  have  been  added  need  to  be  examined  closely. 
This  process  is  given  funher  prionties  by  identification  of  the  important  features  through 
application  of  sensitivity  analysis  and  reaction  flux  analysis  techniques. 

The  validation  at  each  level  of  development  of  a  reaction  mechanism  requires  experimen¬ 
tal  data  from  several  sources  such  as  shock  tubes,  flow  reactors,  and  static  reactors.  It  is 
not  sufficient  to  test  a  mechanism  by  comparison  with  a  single  experiment,  because  differ¬ 
ent  elementary’  reactions  can  be  dominant  under  different  experimental  conditions.  For  ex¬ 
ample.  with  many  hydrocarbons,  reactions  between  fuel  molecules  and  H-atoms  are 


dominant  in  fuel-rich  conditions,  but  much  less  imponant  in  fuel-lean  conditions  where 
reactions  between  fuel  and  O  and  OH  radicals  are  most  important.  Furthermore,  no  single 
kinetics  expenment  covers  the  entire  range  of  pressure/iemperature/ residence  time  pa¬ 
rameters  which  will  be  important  under  combustion  conditions.  Some  reactions  are  negli¬ 
gible  except  under  high  temperature  shock  tube  conditions,  while  others  become 
unimportant  for  temperatures  above  1000  K.  Unless  a  mechanism  is  to  be  used  only  for  a 
restricted  class  of  applications  under  a  specific  set  of  system  parameters,  it  must  be  vali¬ 
dated  by  comparison  with  expenmental  data  over  wide  ranges  of  physical  conditions. 

The  three  experiments  mentioned  above  decouple  the  chemistry  from  transport  phenome¬ 
na  and  excessive  heat  release  rates,  and  provide  mechanistic  kineoc  data  for  different  time 
scales  of  reaction.  Typically,  shock  tubes  cover  1-1000  |i.s  and  hence  are  used  to  study 
high  temperature  kinetics  (T  >  1000  K).  Row  reactors  have  time  scales  typically  between 
10  and  1000  ms  and  are  used  to  study  kinetics  at  intermediate  temperatures  (800  K  <T  < 
1200  K)  while  static  reactors  have  a  time  scale  between  0.1  and  KXX)  s  and  are  used  to 
study  kinetics  at  low  temperatures  (T  <  800  K). 

Tlie  approach  for  mechanism  development  discussed  above  for  hydrocarbon  oxidation  is 
also  applicable  to  the  construction  of  reaction  mechanisms  for  the  gas-phase  intermediates 
of  propellant  combustion.  An  imponant  component  to  this  process  not  discussed  here  is 
the  experimental  and  theoretical  studies  on  individual  species  and  reacnons.  descnbed  in 
the  bottom  half  of  FIGURE  2.  Such  studies  provide  the  elemental  parameter  input  necessary 
for  mechanism  development.  These  two  ty  pes  of  kinetic  studies,  that  is.  those  devoted  to 
complex  reactions  and  those  devoted  to  isolated  reactions,  must  be  conducted  interactively 
for  efficient  and  successful  model  development.  For  example,  the  rates  of  all  possible  ele¬ 
mentary  reactions  cannot  be  practically  measured  in  the  laboratory;  however,  mechanistic 
studies  combined  with  sensitivity  analysis  may  be  used  to  pinpoint  and/or  to  place  priority 
on  those  individual  reactions  w  hich  require  the  most  attention. 

Once  a  gas-phase  kinetics  model  has  been  developed,  validation  of  coupled  gas-phase  heat 
release  rates  and  gas-phase  transpon  is  required  to  complete  the  development  and  valida¬ 
tion  of  the  entire  gas-phase  thermochemical,  kinetics  and  transpon  model.  The  procedure 
is  again  hierarchical  and  is  show  n  schematically  in  FIGURE  3.  Data  from  premixed  flames 
are  necessary  to  validate  heat  release  rates  and  flame  speed  predictions.  In  premixed 
flames,  kinetics  and  transport  are  of  almost  equal  importance;  how  ever,  the  system  is  al¬ 
most  entirely  dnven  by  the  heat  release  and  hence  the  temperature  profile  through  the 
flame.  The  ability  to  deconvolute  the  kinetics,  which  produce  this  heat  release,  is  very  dif¬ 
ficult  due  to  simultaneous  transport  processes  and  the  high  sensitivity  of  measured  observ¬ 
ables  to  the  temperature  measurements.  Data  from  diffusion  flames  are  necessary  to 
validate  transport  and  extinction  phenomena,  in  addition  to  gas-  phase  kinetics  and  heat 
release  rates. 

The  objectives  of  this  research  program  are  to  obtain  kinetic  data  from  flow  reactor  exper¬ 
iments  and  to  develop  and  validate  reaction  mechanisms  for  gas-phase  intermediates  im¬ 
portant  to  the  combustion  of  energetic  polycyclic  nitramines  and  energetic  binders. The 
inherent  nature  of  the  research  requires:  ( 1)  an  experimental  effort  for  determining  reac¬ 
tant.  intermediate,  and  product  species  profiles  as  a  function  of  reaction  time  m  a  flow  re- 
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FIGURE  3.  Hierarchical  Development  and  Validation  ol  Gas  Phase 
Model 


actor.  (2)  comprehensive  reaction  mechanism  development  and  accompanying  computer 
modeling.  (3)  application  of  existing  and  developing  sensitivity  analysis  techniques  to  ex¬ 
tract  the  maximum  information  from  the  first  two  parts  of  the  research.  (4)  a  research  ef- 
fon  coordinated  and  integrated  with  other  laboratories  studying  similar,  yet 
complimentary,  gas-phase  kinetics,  and  (5)  coordinated  efforts  with  other  laboratories  to 
interface  the  gas-phase  mechanisms  developed  with  heterogeneous  reaction  chemistry  and 
solid  propellant  surfaces. 


1.2  Program  Overview  and  Status 


The  experimental  asnecis  of  this  work  are  being  conducted  in  a  variable  temperature  (600 
K<T  <  1200  *0,  vanable  pressure  ( 1  atm  <  P  <  18  attn)  flow  reactor  with  residence  times 
ranging  from  10  to  10  5.  A  second  variable  temperature,  atmosphenc  pressure  reactor 
(800  K  <T  <  1200  K)  is  also  available  for  use.  The  measurement  of  reactant,  intermediate, 
and  product  species  concentrations  as  a  function  of  reaction  time  via  continuous  sampling 
and  thermal  quenching  provide  not  only  substantial  definition  of  the  phenomenology  of 
reaction  mechanisms,  but  a  much  more  constrained  set  of  pure  kinetic  data  than  can  be  de¬ 
rived  from  flame  and  shock  tube  experiments.  Analytical  detection  and  quantification  of 
stable  species  includes  on-line  analysis  by  Fourier  transform  infrared.  In  addition.  CO  and 
CO2  concentrations  are  measured  continuously  by  non-dispersive  infrared.  NO  and  NO2 
concentrations  by  chemiluminescence,  and  O2  concentrations  by  thermomagnetic  analy¬ 
sis.  Off-line  analysis  of  stable  species  includes  gas  chromatography  and  gas  chromatogra¬ 
phy/mass  spectrometry.  Measurement  of  OH  radical  concentrations  are  obtained  by  line  of 
site  absorption  and  laser  induced  fluorescence  techmques.  Temperature  measurements  are 
obtained  by  fine  wire  silica  coated  Pt6%Rh/Pt30%Rh  thermocouples.  Velocity,  and  hence 
relative  reaction  residence  times,  are  obtained  from  reacting  flow  pitot  static  probe  mea¬ 
surements  and  non-reacting  flow  hot  wire  measurements. 

The  modeling  aspects  of  the  program  emphasize  the  use  of  hierarchical  reaction  mecha¬ 
nism  construction  and  validation  along  with  reaction  flux  and  gradient  sensitivity  analysis. 
Sensiuvity  analysis  is  used  to  identify  key  species  and  reactions,  and  hence,  provide  an  in¬ 
teractive  tool  for  coordinating  experiment  and  model  development. 

Kinetic  mecharustic  studies  are  utilized  for  understanding  the  expenmental  observations 
performed  as  part  of  this  program  as  well  as  for  developing  comprehensive  reaction  mech¬ 
anisms  for  use  in  a  general  propellant  combustion  model.  Construction  of  comprehensive 
reaction  mechanisms  is  to  be  carried  out  in  conjunction  with  other  gas-phase  experimental 
and  modeling  programs  which  use  shock  tubes,  premixed  flames,  and  diffusion  flames  as  a 
source  of  kinetic  data.  Each  of  these  experiments  provide  kinetic  data  that  is  complimenta¬ 
ry  to  flow  reactor  data,  which  collectively  provide  the  necessary  data  for  developing  and 
validating  the  most  complete  comprehensive  model,  for  use  in  both  homogeneous  and 
transpon-coupled  situations. 

The  chemical  systems  of  interest  include  the  fuels.  H2.  CO.  CH2O.  CH4.  C2H4.  HCN  and 
NH3  and  oxidizers  NO2.  N2O,  and  NO.  Flow  reactor  studies  will  be  carried  out  with  dilute 
combinations  of  fuel  and  oxidizer  chosen  from  this  list,  reacung  in  a  bath  of  nitrogen.  The 
selection  of  these  fuels  and  oxidizers  is  based  on  existing  preliminary  expenmental  and 
theoretical  evidence  that  HCN.  CH2O,  N2O  and  NO2  are  major  gas-phase  intermediates  of 
polycyclic  nitramine  combustion,  each  being  produced  near  or  at  the  surface  of  the  burn¬ 
ing  propellant.  Ammonia  is  aconvement  source  of  NH2  and  NH  radicals,  both  of  which 
are  important  unstable  intermediates  in  the  oxidation  and  decomposition  of  nitramines. 
Ethylene  and  methane  are  selected  as  representative  gas-phase  intermediates  produced 
from  the  combustion  of  energetic  binders.  The  proposed  hierarchical  structure  among  the 
fuels  follows  three  paths;  a  reaction  sequence  of  stable  carbon  species.  C2H4  ->  CH4  -> 
CH2O  ->  CO/H;0  ->  H2.  and  two  reaction  sequences  of  stable  fuel  bound  nitrogen  spe- 
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cies.  HCN  ->  C0/H20->  H2  and  NH3  ->  H2.  Here  the  kinetics  of  each  system  to  the  right 
of  the  other  is  a  submechanism  of  the  former.  The  oxidizer  mtrogen  chemistry  typically 
follows  the  sequence  NO2  ->  NO  ->  N2O  ->  N2.  Studies  will  begin  with  the  simplest  sys¬ 
tems  and  proceed  towards  the  more  complex. 

The  result  of  this  work  will  be  state-of-the-art  detailed  reaction  mechanisms  validated 
over  the  temperature  and  pressure  ranges  of  the  flow  reactor  experiments  for  the  chemical 
systems  described  above.  With  results  from  other  gas-phase  kinetic  experiments,  a  single 
comprehensive  reaction  mechanism  may  be  developed  for  use  in  1-D  and  2-D  composite 
propellant  combustion  models  which  include  coupled  gas-phase/heterogeneous  chemistry. 

During  the  past  year,  a  Nicolet  Model  730  FTIR  was  installed  and  calibrated  for  the  mea¬ 
surement  of  oxides  of  nitrogen,  simple  hydrocarbons,  carbon  oxides  and  water.  Under  oth¬ 
er  funding,  a  continuous  selective  analytical  technique  has  also  been  developed  for 
molecular  hydrogen.  Experiments  and  modeling  calculations  were  conducted  on  the  kinet¬ 
ics  ofN20  and  A'(92  decomposition  (1000  K  -1100/(.  3  attn.  -tann.).  and  on  the  reactions 
of //2  with  N2O  (925  K  -  1073  A",  1  ann.  -  3  attn.)  and  Hi  with.V02  (803  K  -  873  K.  3  attn.) 

1.3  Report  Outline 

Section  2.0  presents  the  flow  reactor  experiment.  A  discussion  of  the  modifications  made 
to  this  system  that  were  required  in  order  to  study  the  kinetics  of  nitrogen  oxidizers  is  giv¬ 
en.  as  well  as  the  procedures  used  to  measure  and  quantify  the  concentrations  of  the  oxides 
of  mtrogen.  A  discussion  of  the  modeling  approach  is  given  in  Section  3.0.  The  current  re¬ 
action  mechanism  is  also  described  in  this  section.  In  Section  4.0.  a  description  of  the 
studies  on  N2O  decomposition  and  on  its  reaction  with  H2  are  given.  Section  5.0  presents 
the  results  of  kinetic  studies  on  NO2  decomposition  and  its  reaction  with  H2.  Finally.  Sec¬ 
tion  6  describes  the  current  and  on-going  research. 

2.0  FLOW  REACTOR  EXPERIMENT 

2.1  The  Atmospheric  and  Variable  Pressure  Flow  Reactors 

Flow  reactor  experiments  are  unique  in  providing  experimental  data  in  a  range  of  condi¬ 
tions  typically  not  accessible  in  shock  tubes  and  staac  reactors  yet  still  of  significant  im¬ 
portance  to  propellant  combustion.  From  such  expenments.  the  evolution  of  stable 
reactants,  intermediates  and  product  concentrations  for  chemical  reacting  systems  are  de¬ 
fined  in  considerable  detail  producing  highly  constrained  sets  of  data  for  model  develop¬ 
ment.  In  fact,  this  technique  (described  in  detail  in  APPENDIX  A)  has  been  the  source  of 
much  of  the  previous  intermediate  temperature,  atmosphenc  pressure  data  for  the  oxida¬ 
tion  of  methane,  ethylene,  propane,  n-butane.  carbon  monoxide/hydrogen,  octane  (Haut- 
man.  1981:  Proscia.  1983;  Westbrook  and  Pitz.  1984;  Pitz,  etal..  1984;  Westbrook  et  al.. 
1983;  Westbrook  and  Dryer,  1984;  Yetter.  1985;  Brezinsky  and  Dryer.  1985,  1986;  Dryer 
and  Brezinsky,  1986;  Axelsson  et  al.,  1988;  Dryer.  1989)  and  simple  aromatics  (Brezin¬ 
sky.  198.*';  and  subsequent  work  of  Prof  1  Glassman  and  his  colleagues').  These  data  have 


provided  not  only  phenomenological  insights  into  mechamsms  for  pyrolvsis  and  oxida¬ 
tion,  but  a  source  of  information  for  validation  of  detailed  mechanisms,  and  a  means  to 
study  the  intimate  coupling  of  elementary  reactions. 

The  variable  pressure  flow  reactor  accesses  a  range  of  pressures  (0. 1  to  28  atmospheres  by 
design,  1.0  to  18  atm.ospheres  in  present  operation),  and  a  wider  range  of  temperatures 
(from  as  low  as  600  K  to  1200  K)  with  improved  experimental  control  and  diagnostics,  in 
companson  to  those  of  the  atmosphenc  pressure  flow  reactor  (appendix  A).  Available  re¬ 
action  times  as  shon  as  50  milliseconds  and  as  long  as  10.000  milliseconds  are  possible, 
all  with  the  same  experimental  method,  and  without  having  to  consider  heterogeneous 
wall  effects  on  the  measured  homogeneous  kinetics.  In  addition  to  the  pressure  range  now 
accessible  in  this  system,  several  other  new  expenmental  features  have  been  incorporated 
into  this  system.  The  addition  of  a  Nicolet  Model  730  FTIR  permits  on-line  measurement 
of  most  of  the  IR  absorptive  stable  species  of  interest  under  this  research  effort  to  deter¬ 
mine  levels  of  about  0. 1  pp>»-  The  new  design  of  the  high  pressure  flow  reactor  also  allows 
for  improved  optical  access,  permitting  simultaneous  insitu  optical  diagnostic  detection  of 
unstable  species.  The  integration  of  the  entire  experiment  (experimental  operation  and 
data  acquisition)  into  a  closed  loop  control  computer  based  system  allows  for  more  fre¬ 
quent  expenments  and  improved  data  collection. 

The  upper  residence  time  limit  cntena  for  flow  reactors  is  the  maintenance  of  convective 
residence  times  which  are  short  in  comparison  to  diffusion  times  to  the  reactor  wall  sur¬ 
face.  This  limiting  criteria  is  strongly  influenced  by  the  selecnon  of  the  flow  reactor  diam¬ 
eter.  Estimating  the  order  of  the  diffusion  time  to  the  wall  as 

Td  =  d^  /4  D 

where  d  is  the  diameter  of  the  reactor  tube  and  D  is  the  mass  diffusion  coefficient,  and  the 
convective  residence  time  as 

\  =  L/u 

where  L  is  the  length  of  the  reactor  duct  and  u  is  the  mean  flow  velocity,  »  T;.  when 

iRe)(Sc) »  4L/d 

where  Re  designates  the  flow  Reynolds  number  based  on  the  reactor  duct  diameter  and  Sc 
designates  the  Schmidt  number  of  the  carrier  gas.  For  the  present  parameter  conditions, 
this  translates  to 

(Re)(Sc) »  40 

As  the  Schmidt  number  is  the  order  of  one.  this  means  the  flow  Reynolds  numbet  should 
be  much  larger  than  40  to  meet  this  condition.  The  reactor  is  intended  to  be  run  only  at 
Reynolds  numbers  greater  than  about  1000.  As  the  duct  L/d  downstream  of  the  mixer  sec¬ 
tion  will  always  be  less  than  10.  the  flow  is  not  fully  developed.  As  the  boundary  layer  de¬ 
velops  at  the  duct  surface,  radially  uniform  conditions  continue  to  exist  near  the  flow 
centerline,  .As  the  flow  field  further  develops,  axial  velocity  will  change  as  a  function  of 


distance  from  the  mixer  section.  Axial  velocity  profiles  must  therefore  be  defined  in  order 
to  determine  the  relative  reaction  time  as  a  function  of  distance  from  the  mixer  to  the  sam¬ 
pling  position. 

However,  very  important  implications  result  from  these  observations,  particularly  at  in¬ 
creased  operating  pressures.  FIGURES  4  and  5  show  Reynolds  number  of  the  flow  for  a  10 
cm  diameter  reactor  duct  as  a  function  of  available  reaction  ume  for  various  operating 
pressures  and  temperatures.  Note  that  at  higher  pressures  and  lower  reaction  temperatures, 
the  reaction  test  times  at  Reynolds  numbers  exceeding  1000  cover  a  range  from  10  milli¬ 
seconds  to  more  than  10.000  milliseconds.  This  new  facility  therefore  presents  the  possi¬ 
bility  of  studying  kinetic  issues  over  a  very  range  of  conditions,  including  a  temperature 
range  previously  available  only  in  static  reactors. 

One  of  the  basic  design  pnnciples  of  static  reactors  is  that  diffusive  transpon  is  rapid  with 
regard  to  chemical  reaction  and  energy  release  so  that  the  local  conditions  within  the  reac¬ 
tor  remain  homogenous,  and  isothermal.  Reaction  times  under  these  conditions  are  typi¬ 
cally  in  excess  of  many  tens  of  seconds,  and  thus,  diffusion  of  radical  and  unstable 
molecular  intermediates  to  surfaces  result  in  perturbations  of  the  gas  phase  chemistry. 
There  are  no  such  surface  effects  in  the  variable  pressure  flow  reactor. 

Tlie  design  concept  of  the  variable  pressure  flow  reactor  is  based  on  fixing  the  diagnostic 
sampling  position,  and  moving  the  point  of  fuel  injection  relative  to  this  location  to  vary 
the  relative  reaction  time  one  wishes  to  observe.  The  approach  is  very  similar  to  that 
which  is  normally  performed  in  the  study  of  pre-mixed  laminar  one  dimensional  flames 
where  the  burner  rather  than  the  diagnostic  sampling  position  is  moved.  This  approach  not 
only  permits  very  shon  gas  sampling  residence  times  to  continuous,  on-line  diagnostic  in¬ 
struments  (important  for  on-line  measurement  of  low-stability  molecular  species  such  as 
aldehydes  and  other  oxygenates),  but  simple  accommodation  of  optical  diagnostics  at  the 
sampling  location.  Crossed-beam  optical  access  ports  (normal  to  the  flow  direction)  are 
positioned  at  the  same  location  at  w  hich  a  hot  water  cooled,  wall  convection  quenched, 
gas  sampling  probe  and  silica  coated  thermocouple  probe  are  axially  located.  (Optical  diag¬ 
nostic  measurements  of  OH  using  exisung  line  resonance  absorption  or  fluorescence  tech¬ 
niques  (Linteris.  1991)  can  be  much  more  easily  accomplished  than  in  the  aimosphenc 
pressure  flow  reactor,  and  other  insitu  optical  measurements,  including  insitu  FTIR  and 
other  resonance  absorption  measurements  can  also  be  added  to  the  facility  in  the  future. 

A  schematic  diagram  of  the  variable  pressure  flow  reactor  is  given  in  FIGURE  6.  Nitrogen 
earner  gas  with  oxidant  flows  from  left  to  right  in  the  figure.  Fuel,  along  with  about  10%of 
the  total  mtrogen  carrier  flow  is  injected  into  the  carrier  and  oxidizer  stream  at  a  mixer  on 
the  left  end  of  a  conical  silica  foam  diffuser.  The  mixer  is  jet  stirred  and  has  a  nominal 
turnover  time  of  about  0.59c  of  the  total  test  time.  The  flow  reactor  is  operated  a'  a  steady, 
isobanc  flow  device  and  the  ume  from  mixing  to  sample  position  is  vaned  by  moving  the 
mixer  /diffuser  assembly  along  the  quanz  test  tube.  At  the  right  end  is  the  fixed  sample  lo¬ 
cation. 

The  shell  is  made  of  1 2"  standard  weight  carbon  steel  pipe  and  is  ASME  stamped  for  op¬ 
eration  from  full  vacuum  to  30  arm  and  245  K  to  533  K.  The  pressure  shell  is  designed  to 


FIGURE  4.  Operational  Reaction  Test  time  as  a  Function  of  Flow 
Reynolds  Number  and  Operating  Pressure  for  a  10  cm  Diameter 
Reactor  Duct 


enclose  not  only  the  reactor  section  (the  reactor  duct  heaters  and  insulation),  but  also  the 
positioning  mechanism  and  fuel  vapor  iniection  probe  which  axially  locates  the  fuel/cam- 
er  mixing  section  relative  to  the  sampling  position.  This  design  results  in  very  low  thrust 
to  move  probes  (due  to  pressure  drop  across  flanges)  and  eliminates  the  need  for  critical 
dynamic  seals.  Since  all  flanges  are  at  lower  temperature,  there  is  also  no  need  for  exotic 
materials  of  construction  to  be  used  for  the  pressure  shell  nor  elaborate  designs  required  to 
maintain  pressure  seals. 


The  reactor  duct  itself  is  a  173  cm  long.  2  mm  wall  thickness.  10. 16  cm  inside  diameter, 
fused  silica  tube  which  extends  from  a  Inconel  mating  flange  at  the  fuel  injection  probe 
entrance  to  a  similar  mating  flange  located  at  the  entrance  to  the  exhaust  port  of  the  reactor 


log(t/s) 


FIGURE  5.  Operational  Reaction  Tune  as  a  Function  of  Flow 
Reynolds  Number  and  Flow  Temperatures  and  15  Atmosphere 
Operating  Pressure  for  a  10  cm  Diameter  Reactor  Tube 


shell.  A  maximum  of  about  100  cm  and  a  minimum  of  1  cm  of  axial  distance  from  the  exit 
of  the  mixing  section  to  the  sampling  position  are  possible.  In  principle,  reactor  duct  diam¬ 
eter  should  be  able  to  be  reduced  to  as  small  as  5  cm  inside  diameter  by  replacing  the  in¬ 
conel  mating  flanges,  the  reactor  duct  itself,  and  installing  a  smaller  diameter  mixing 
section.  The  carbon  steel  shell  is  manufactured  in  two  separate  parts  in  order  to  accommo¬ 
date  such  changes. 

A  cross  secuon  of  the  fuel  mixer  and  diffuser  is  shown  in  FIGURE  7.  The  mixer  section  is 
approximately  47  cm  long  and  is  machined  from  a  low  porosity,  silica  foam  block  to  an 
outside  diameter  about  0. 1  mm  less  than  the  nominal  10.16  cm  inside  diameter  of  cylindri¬ 
cal  reactor  duct.  The  interior  contour  of  the  mixer  section  is  shaped  as  a  converging/di- 


Cross  Section  of  VPFR  Vessel  and  Internals  (to  scale) 


FIGL'RE  6.  Schematic  of  variable  pressure  flow  reactor  design 


verging  nozzle  with  a  throat  diameter  of  2.5  cm  and  an  expansion  half  angle  of  7  degrees. 
Flow  resistance  through  the  long  annular  space  between  the  mixer  section  and  reactor  duct 
wall  is  very  large  in  comparison  to  that  through  the  central  venturi  and  nearly  all  of  the 
earner  flow  therefore  passes  through  the  central  venturi.  One  inch  long  ceramic  fiber  pis¬ 
ton  ring  seals  are  added  at  both  the  upstream  and  downstream  ends  to  mimmize  leakage 
between  the  diffuser  and  the  quartz  test  tube.  The  mixer  -  diffuser  section  consists  of  a 
central  round  baffle  plate  which  forces  the  nitrogen-oxidizer  flow  out  to  the  reactor  tube 
wall  and  then  radially  inward  through  a  gap  between  the  plate  and  the  upstream  end  of  the 
diffuser.  Fuel  and  diluent  mtrogen  are  delivered  to  the  mixer  through  a  fuel  vapor  injection 
probe,  which  is  used  also  to  move  the  diffuser.  The  fuel  vapor  injection  probe  is  partially 
insulated  from  the  carrier  flow  to  prevent  premature  pyrolysis  reactions  of  the  niuogen  di- 


FIGURE  7.  Schematic  of  mixer  and  diffuser  section 


luted  fuel  vapor  inside  the  probe  from  occurring.  All  fuels  are  premixed  into  varying 
amounts  of  added  nitrogen  diluent  to  reduce  residence  time  within  the  fuel  injection  probe. 
The  injector  for  the  present  expenments  was  fabricated  from  a  3 1 6  stainless  steel  6.35  mm 
outside  diameter  X  0.9  mm  thick  wall  tube.  The  fuel-diluent  flow  is  then  injected  radially 
outw  ard  into  this  radially  inw  ard  flow  through  eight  0.25  mm  diameter  jets  equally  spaced 
around  the  circumference  of  the  injector  tube.  All  gas  flows  are  measured  with  Teledyne- 
Hastings  mass  flow  meters  with  a  maximum  range  of  uncertainty  in  any  gas  flow  estimat¬ 
ed  to  be  about  2%  of  full  scale,  which  for  the  present  experiments  translates  into  an  uncer¬ 
tainty  range  of  not  more  than  b'T  for  the  nitrogen  and  oxidizer  flows. 


FIGURE  8.  Block  Diagram  of  the  Control  System  for  the  Variable 
Pressure  Flow  Reactor 


Because  of  the  relatively  high  heating  rates  involved  (tens  of  h\-),  as  well  as  the  tempera¬ 
tures  involved  in  some  of  the  systems  to  be  studied  (as  high  as  1.100  C.  the  limitation  of 
the  quartz  test  tube),  an  arc  plasma  gas  heater  is  used  to  heat  the  carrier  gas.  The  arc  plas¬ 
ma  torch  heater  is  similar  to  that  used  for  the  last  fifteen  years  at  Pnnceion.  A  new  micro¬ 
processor-based  flow  and  arc  power  control  system  has  been  designed  to  vary  the  plasma 
power  input  and  thermal  efficiency  for  controlling  the  earner  flow  temperature,  as  moni¬ 
tored  at  the  fuel/carrier  mixer  location  (thus  radiation  losses  to  the  fuel  injection  probe 
from  the  carrier  gas  do  not  affect  the  inlet  earner  temperature).  Automated  control  systems 
represent  another  imponant  design  improvement,  since  the  precision  of  the  quasi-steady 
operation  of  a  flow  reactor,  parucularly  its  temperature  field,  affect  the  accuracy  of  the  ki¬ 
netic  measurements.  To  our  knowledge  other  flow  reactor  systems  have  not  yet  developed 


FIGLRE  9.  Gas  Flow  Control  Schematic  for  the  Variable  Pressure 
Flow  Reactor 


or  implemented  integrated  controls.  A  35  kw  electric  resistance  gas  heater  has  been  devel¬ 
oped  to  replace  the  arc  plasma  gas  heater  to  operate  at  low  flow  (less  than  5  g/s  total),  low 
duty  (less  than  10  high  pressure  (above  10  arm)  conditions.  ncuRE  8  shows  schemat¬ 
ically  the  control  approach  used  in  the  new  reactor  facility  design,  w'hile  FIGURE  9  shows 
the  flow  systems  design.  All  flows  are  monitored  by  electronic  mass  flow  meters  and  posi¬ 
tioning  of  motor  driven  metering  valves  are  controlled  by  Proportional  -  Integral  -  Deriva¬ 
tive  (PID)  Controllers  which  sense  the  mass  flow  meter  output.  Similar  closed  loop 
controllers  are  used  to  control  all  heater  temperatures  and  the  flow  reactor  operating  pres¬ 
sure.  Operating  pressure  of  the  reactor  is  controlled  by  closed  loop  control  of  a  high  tem¬ 
perature  back  pressure  valve  in  the  exhaust  flow  secdon  The  PID  controllers  may  be 
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nCL'RE  10.  Normalized  centerline  velocity  profile  from  hot  wire 
anemometry  measurements 


manually  operated.  Each  controller  is  also  linked  to  a  IBM/PC/AT  microprocessor  via  a 
RS-485  pon  connection  so  that  all  control  parameters  may  be  read  or  altered  by  the  PC. 

In  order  to  translate  flow  reactor  position  into  reaction  time,  the  mean  velocity  distribution 
along  the  centerline  is  required.  The  centerline  velocity  profile  is  measured  with  a  TSl 
Model  1"50  hot  w  ire  anemometer.  These  measurements  are  performed  in  the  reactor  un¬ 
der  cold  flow  conditions.  Typical  mean  velocity  measurements  are  reponed  in  FIGURE  10. 
The  corresponding  root  mean  square  velocities  reported  in  FIGURE  1 1 .  The 

U  rms  measurements  show  that  the  turbulence  is  decaying  over  the  length  of  the  reactor. 
However,  examination  of  FIGURE  10  shows  that  the  flow  in  the  diffuser  and  test  section  is 
nearly  plug  flow.  The  reaction  time  is  obtained  by  integration  of  the  test  section  velocity 
profile.  An  effective  flow  area  table  for  the  diffuser  and  test  section  is  constructed  from  the 
cold  flow  hot  wire  anemometry  velocity  survey.  Effective  flow  areas  from  this  table,  along 
with  the  measured  flow  rates,  temperatures  and  pressures  are  used  to  obtain  the  local  val¬ 
ues  of  u(x).  Residence  time  is  then  obtained  by  numerically  integraung 
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The  maximum  uncenainty  in  reaction  ume.  including  the  mixing  time  (0.5%).  the  uncer¬ 
tainty  in  the  velocity  survey  (6%),  the  uncenainty  in  flow  measurements  (6%),  the  uncer¬ 
tainty  in  probe  position  (0.4%),  and  the  sample  quench  time  is  +/-  7%. 


The  sample  temperature  is  measured  w  ith  a  silica  coated  Pt6%Rh/Pt30%Rh  thermocou¬ 
ple.  The  total  estimated  uncenainty  in  temperature  measurement,  including  allowance  for 


FIGURE  11.  Centerline  turbulent  intensity  profile  from  hot  wire 
anemometry  measurements 


radiation  errors  due  to  the  proximity  of  the  cold  sample  probe  is  estimated  at  +/■  lO/w.  The 
relative  uncertainty  in  temperatures  within  an  experiment  are  estimated  to  be  +/-  1  K. 

A  hot  water  cooled,  stainless  steel  sample  probe  is  mounted  in  the  reactor  end  flange  and  is 
used  to  continuously  extract  and  convective  quench  a  small  portion  of  the  reacting  gases. 
Quench  times  estimated  for  atmospheric  pressure  flow  reactor  operation  at  1  atm  are  of  the 
order  of  0.0005  s  for  quench  to  600  K  to  0.002  s  for  quench  to  400  K  (not  more  than  2%  of 
total  test  time).  These  estimates  are  based  on  a  finite  difference  model  that  was  calibrated 
with  actual  temperature  measurements  (Emdee.  1988).  The  quench  time  scales  with  the 
characteristic  time  for  diffusion  of  energy  from  the  centerline  of  the  sample  tube  to  the 
wall  of  the  sample  tube,  and  thus,  the  quench  time  is  proponional  to  the  experimental 
pressure.  Because  kinetic  times  are  typically  inversely  proportional  to  the  pressure,  the 
quench  times  are  estimated  to  be  of  the  order  of  1%  of  the  test  times  at  all  pressures.  From 
the  probe,  the  gas  sample  flow  is  sent  to  a  FTIR.  a  continuous  electrochemical  analyzer  for 
O2.  continuous  NDIR  analyzers  for  CO  and  CO2,  and  to  a  continuous  selective  detector 
for  H2.  A  diagram  of  the  sample  tram  is  given  in  rgure  12.  Data  acquired  by  these  in¬ 
struments  IS  forwarded  to  an  a/d  board  interface  and  recorded  on  a  hard  disk  based  micro¬ 
processor  system.  As  in  the  atmospheric  pressure  flow  reactor,  this  continuous  sampling 
and  analysis  system  will  also  accommodate  intermittent  sample  trapping  for  gas  chro¬ 
matographic  analysis.  A  gas  sampling  and  storage  technique  for  interfacing  with  gas  chro¬ 
matographs  and  gas  chromatograph/mass  spectrometers,  developed  here  for  the  AFPR  and 
used  extensively  in  our  previous  work,  continues  to  be  available. 


Sample/Calibration  Gases 


CO2  NDIR 
Horiba 


Exhaust 


r 


_! _ 

j 

H2  Meter 


Exhaust 


FIGL’RE  12.  Schematic  of  sample  train 


Species  and  temperature  profiles  are  obtained  by  positioning  the  diffuser  -  mixer  section  at 
equally  spaced  locations  along  the  length  of  the  reactor  centerline.  About  2.5  minutes  of 
time  IS  required  for  the  sample  flow  to  purge  the  0.5  liter  FTIR  gas  cell.  Stabilization  of 
the  oxygen,  carbon  monoxide,  carbon  dioxide  concentrations  as  indicated  by  the  or.  line 
analyzers,  which  are  upstream  of  the  FTIR  in  the  sample  hne.  provides  venfication  of  ade¬ 
quate  purging.  The  FTIR  purge  time  is  currently  the  factor  that  limits  the  data  collection 
rate. 


The  vanable  pressure  flow  reactor  facility  is  located  immediately  adjacent  to  a  laser  diag¬ 
nostic  facility  (developed  under  a  DOE  University  Research  Initiative  grant)  w  hich  has  al¬ 
ready  been  used  to  make  laser  resonance  absorption  and  fluorescence  measurements  of 


OH  in  the  atmospheric  pressure  reactor  (Lintens.  1991).  Adaptation  of  this  facility  to  the 
vanable  pressure  flow  reactor  involves  primarily  optics  development  for  uansfening  the 
laser  beam  to  the  facility  and  mating  window  designs  for  the  reactor  shell,  and  should  be 
completed  within  the  coming  year. 

2.2  Flow  reactor  modifications  for  N2O  and  NO2 

A  number  of  modifications  to  the  appararas  and  experimental  procedure  have  been  neces¬ 
sary  for  studies  with  N2O  and  NO2  as  oxidizers.  First,  both  N2O  and  NO2  decompose 
slowly  at  the  conditions  under  study  here  (see  the  discussion  later  in  this  report).  Conse¬ 
quently.  their  addition  to  the  carrier  gas  upstream  of  the  mixer-diffuser  section  was  not  ac¬ 
ceptable.  This  difficultly  for  the  present  studies  was  overcome  by  introducing  the  N2O  and 
NO2  through  the  fuel  injector  right  at  the  mixer.  Thus,  when  the  H2  reaction  was  studied, 
H2  w'as  added  upstream  of  the  mixer-diffuser  section. 

Second,  because  of  the  large  mass  flows  of  N2  required  in  the  system  (e.g.  35  g/s).  nitro¬ 
gen  IS  stored  as  a  liquid  and  vaporized  prior  to  its  introduction  into  the  reactor.  Small 
quantities  of  molecular  oxygen  (typically  less  than  20  ppm)  may  be  present  in  the  flow, 
which  will  be  shown  later  to  affect  the  kinetics  of  the  N2O  system.  This  problem  was  re¬ 
duced  by  the  addition  of  the  H2  upstream  of  the  mixer,  allowing  time  for  the  O2  to  be  con¬ 
sumed  by  the  H2  prior  to  N2O  or  NO2  injection. 

The  third  and  most  senous  problem  was  the  present  choice  of  construction  material  uti¬ 
lized  for  the  baffle  plate  and  the  injector  tube  itself  stainless  steel).  Hot  stainless  steel  sur¬ 
faces  can  act  as  a  catalyst  for  NO2  particularly  at  surface  temperatures  higher  than 
approximately  930  K.  For  the  present  set  of  experiments  with  NO2.  this  temperature  was 
an  upper  limit.  The  mixer-diffuser  section  has  been  redesigned  and  all  of  the  stainless  steel 
components  are  currently  being  replaced  with  quartz. 

2.3  Measurement  of  Oxides  of  Nitrogen  and  Water  Vapor 

A  Nicolet  Model  730  founer  transform  infrared  analyzer  (FTIR)  is  used  to  measure  the 
oxides  of  nitrogen  and  water  vapor  for  the  current  results  and  will  be  used  for  measure¬ 
ment  of  other  reactants  in  future  studies,  e.g.,  compounds  containing  up  to  four  carbon  at¬ 
oms.  On-line  FTIR  should  also  be  able  to  measure  accurately  some  species  which  cannot 
be  stored  (without  chemical  stabilization)  for  later  off-line  analysis  using  gas  chromato¬ 
graphs.  etc.,  such  as  hydrogen  peroxide,  formaldehyde,  and  other  oxygenates. 

The  analyzer  has  a  0.5  liter  gas  cell  with  a  optical  path  length  of  8  m.  A  liquid  nitrogen- 
cooled  MCTA  detector  is  used  to  measure  spectral  absorbance  with  a  resolution  of  0.5  cm' 
'.  Multivariate  least  squares  fit  software  is  used  for  quantification  of  the  spectral  data. 

Calibrations  for  NO.  N2O  and  NO2  were  obtained  by  measuring  absorbance  for  different 
concentranons  of  each  compound  in  nitrogen  (minimum  of  three)  from  calibrated  gas  cyl¬ 
inder  bottles  supplied  by  MATHESON.  Examples  of  spectral  calibration  data  are  given  in 
FIGURES  13-15  for  7260 ppm  of  N20in  N2. 500  ppm  NO  in  N2.  and  2996 ppm  of  NO2  in 
N;.  respectively. 
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FIGURE  13.  Absorbance  Spectrum  from  Nitrous  Oxide  in  Nitrogen 


The  FTIR  has  also  allowed  measurement  of  water  vapor  concentration.  This  measurement 
IS  important  because  it  permits  closure  of  the  hydrogen  atom  mass  balance,  which  is  typi¬ 
cally  not  available  in  most  kinetic  studies.  Water  calibrations  were  made  from  the  catalytic 
oxidation  of  a  known  amount  of  H2  in  an  O2/N2  mixture.  An  example  of  the  spectra  for 
4000 ppm  of  H2O  in  N2  is  given  in  FIGURE  16. 

An  example  of  the  spectra  from  an  actual  H2/N2O  experiment  is  reponed  in  FIGURE  17 
illustrating  the  well  resolved  spectra.  All  of  the  spectra  reported  here  were  obtained  with 
32  scans. 

2,4  .Measurement  of  Molecular  Hydrogen 

Molecular  hydrogen  is  present  in  most  combustion  processes  in  concentrations  from  pans 
per  million  {ppm)  to  several  percent  by  volume.  However,  hydrogen  is  inherently  difficult 
to  quantify  in  multi-component  mixtures.  The  significance  of  the  hydrogen  molecule,  as 
either  a  pnmary  fuel  or  an  intermediate  in  hydrocarbon  oxidation,  makes  its  measurement 
very  important.  Existing  direct  methods  of  hydrogen  detection  include  gas  chromatogra¬ 
phy  with  a  thermal  conductivity  detector  (TCD);  however  they  require  collection  of  indi¬ 
vidual  gas  samples  and  relatively  long  analysis  times.  Under  other  support,  a  new 
technique  has  been  developed  to  continuously  measure  the  amount  of  hydrogen  in  the 
flow  of  a  gaseous  mixture.  A  detailed  description  of  the  system  will  be  availat’  *.  following 
a  patent  evaluation  on  the  device  (which  is  presently  underway). 
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3.0  MODEL  DEVELOPMENT  AND  VALIDATION 


3.1  Flow  Reactor  Model  and  Sensitivity  Analysis  Equations 

3.1.1  Model  Equations 

In  the  case  of  the  flow  reactor,  convective  transport  is  large  relative  to  the  diffusive  trans¬ 
port  as  a  result  of  operation  with  high  through-put  velocities,  and  thus  the  longitudinal  gra¬ 
dients  are  small.  This  assumption  reduces  the  steady,  one-dimensional  conservation 
equations  of  mass,  species,  and  energy  to  the  following 


M  =  puA 
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In  these  equations  x  denotes  the  independent  spatial  coordinate;  M.  the  mass  flow  rate;  T. 
the  temperature;  the  mass  fraction  of  the  ktfi  species;  p.  the  mass  density;  the  mo¬ 
lecular  weight  of  the  kth  species;  Cp  the  constant  pressure  heat  capacity  of  the  mixture; 
the  specific  enthalpy  of  the  kih  species;  and  A.  the  cross  sectional  area  of  the  reactor  nibe. 


Assuming  M  and  A  as  constants  and  equaung  the  convective  derivative  ud/dx  to  the  un¬ 
steady  term  d  dt.  the  above  equations  can  be  rew  rmen  as. 
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If  the  flow  reactor  is  operated  w  ith  dilute  mixtures,  the  exothermicity  of  the  mixture  is 
small  and  the  system  remains  nearly  isothermal.  Under  these  conditions,  the  energy  equa¬ 
tion  is  eliminated  from  the  above  set. 


3.1.2  Sensitivity  Analysis 

The  role  of  sensitivity  analysis  is  to  probe  the  relationships  between  the  model  parameters 
and  the  solutions  to  the  governing  equations;  that  is.  the  spatial  or  time  dependent  species 
concentrauons,  temperature,  velocity,  etc.  In  essence,  a  sensiuvity  analysis  to  the  above 
flow  reactor  mcxlel  allows  questions  to  be  answered  concerning  (1)  interpretation  and 


analysis  of  models  (2)  relationships  between  model  parameters  and  experiment  and  (3)  in¬ 
version  problems.  Typical  of  these  questions  may  be  the  following: 

•  What  steps  are  sigruficant  in  the  mechanism? 

•  What  role  might  a  particular  observation  play  regarding  identification  or  determination 
of  a  particular  rate  constant? 

•  Can  alterations  in  one  pan  of  a  mechanism  be  compensated  for  by  changes  elsewhere 
while  leaving  the  observations  fixed? 

•  Are  different  observations  in  the  same  chemical  system  truly  independent? 

•  What  physical  content  resides  in  identifiable  features  of  concentration  or  thermal  pro¬ 
files? 

•  How  can  one  establish  the  role  of  initial  species  or  possibly  the  effect  of  the  induction 
of  species  during  the  reaction? 

•  By  what  pathway  through  the  mechanism  does  one  species  control  the  behavior  of  oth¬ 
ers? 

•  How  do  the  various  rate  constants  work  in  concen  to  have  major  or  minor  influences  on 
combustion  phenomena? 

•  Can  the  role  of  missing  reactions  or  steps  in  a  mechanism  be  assessed? 

•  How  may  one  quantitatively  reduce  the  number  of  rate  constants  or  transpon  coeffi¬ 
cients  m  a  reactive  system'’ 

•  Can  systematic  means  be  found  to  lump  combustion  systems'’ 

•  How  do  statistical  errors  in  the  model  input  propagate  to  laboratory  observations? 

•  When  inverting  laborator>'  data  by  residual  minimization,  what  correlabons  are  pro¬ 
duced  among  the  data  and/or  parameters'’ 

The  governing  flow  reactor  equations  described  in  the  last  section  can  be  generalized  as 

F(O.a)  =  0  tEQ6) 

w  here  O  represents  the  vector  of  ,v  dependent  variables,  and  the  vector  a  of  length  M  rep¬ 
resents  the  system  parameters  such  as  activation  energies,  preexponential  factors,  and  oth¬ 
er  quantities  that  enter  the  differential  equations. 

In  a  sensiuvity  analysis,  the  quantities  of  natural  interest  are  the  first  order  sensitivity  coef¬ 
ficients 


do, 

which  provide  a  direct  measure  of  how  the  jih  parameter  controls  the  behavior  of  the  ith 
dependent  variable  at  point  .x  and  time  t.  The  appropriate  equations  for  these  quantities  can 
be  derived  by  differentiating  EQ  4  with  respect  to  .  Noting  that  the  functional  vector  F 


depends  on  a  both  explicitly  and  implicitly  through  the  solution  vector  O,  differentiation 
of  EQ  6  yields 
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{J=]  J . :Vf).  Recalling  that  the  Jacobian  matnx  is  given  by  J  =  BFIBO.  we  have 
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(J=l . M).  Solution  of  these  equations  is  generally  obtained  w  ith  the  simultaneous  solu¬ 

tion  of  the  governing  equation.  In  the  present  work,  the  kinetic  equations  and  sensitivity 
equations  were  solved  using  SENKIN  (Kee  et  ai.  1990) 


The  gradients  studied  in  the  present  work  include  the  sensitivities  of  the  temperature  and 
species  mass  fracuons  with  respect  to  variations  in  the  kinetic  parameters  of  the  mecha- 
msm.  For  example,  the  responses  of  observable  O,  to  variations  in  the  forxs  ard  and  back¬ 
ward  rate  constants.  Ay  and  Ay  are. 

BO,  BO, 

Blnk^  5//iAy 

Because  of  the  principle  of  detailed  balancing,  the  forward  and  reverse  rate  constants  are 
related  through  the  equilibrium  constant.  In  a  sensitivity  analysis,  if  Ay  is  perturbed  and  Ay 
IS  held  constant,  the  equilibnum  constant.  K.  is  simultaneously  varied  by  the  same 
amount.  Posing  the  original  problem  in  terms  of  K  and  Ay.  the  latter  w  hich  we  define  here 
as  A^f,  to  differentiate  it  from  the  onginal  problem  with  Ay  and  Ay  as  the  variables,  yields 
another  set  of  sensitivity  gradients.  These  gradients  are  obtainable  from  the  onginal  set  of 
sensitivities.  It  can  easily  be  shown  that  the  sensitivity  of  observable  (7,  to  A, equals. 

Bo,  Bo,  Bo, 

Blnk„^,  Blnk^  Blrikf^ 

and  the  sensitivity  of  O,  to  the  equilibrium  constant  K  is 

BO,  BO, 

BJnK  ^  ~Blnk^ 

Note  that  if  the  reverse  rate  constant  is  as  important  as  the  forward  rate  constant  and  oppo¬ 
site  in  sign,  the  net  sensitivity  equals  zero.  In  many  situations,  this  is  a  good  indication  that 
the  reaction  is  microscopically  balanced  and  that  the  important  parameter  is  the  equilibn¬ 
um  constant. 


(EQ  10) 


(EQ  11) 


3.2  Chemical  Svsfems  and  the  Modeling  .Approach 

The  proposed  hierarchical  structure  of  chemical  systems  of  imponance  to  propellant  com¬ 
bustion  IS  shown  in  FIGURE  18,  Both  HCN  and  CH2O  are  expected  to  be  major  gas-phase 


intermediates  of  polycyclic  nitramine  combustion.  Decomposiuon  of  RDX  through  con¬ 
certed  detrimenzation  or  imtial  N-NO2  bond  fission  has  been  shown  to  yield  N2O  and 
CH2O  at  low  temperatures  and  HCN  and  NO2  at  high  temperatures.  In  recent  T-jump  ex¬ 
periments  of  Brill  (1991)  for  HMX  decomposition.  \20  and  NO2  were  observed  to  pro¬ 
ceed  all  other  gas  phase  species  from  the  surface.  The  ratio  of  N2O  to  NO2  was  found  to 
decrease  with  increasing  temperatures,  being  approximately  equal  for  temperatures  of 
300-350°C.  which  is  near  the  surface  temperature  of  a  burning  propellant.  Although  the 
ratio  of  CH2O  to  HCN  was  not  observed  to  follow  this  trend,  both  HCN  and  CH2O  are  im¬ 
portant  gas-phase  reactants.  They  are  also  important  intermediates  in  RDX  combustion  as 
shown  in  ncURE  19  for  the  0.5  atiti  pressure  RDX  flame  model  results  of  Melius  (1985). 
Note  from  FIGURE  19,  that  CH2O  is  consumed  very  rapidly  and  prior  to  the  gas-phase 
temperature  rising  above  1200  K  (not  shown  in  this  figure).  Furthermore,  note  that  both 
CH2O  and  HCN  are  nearly  consumed  w'lthin  1  mm  from  the  flame  surface.  Both  of  these 
facts  emphasize  the  importance  of  flow  reactor  experiments,  i.e..  the  chemistry  near  the 
propellant  surface  is  well  within  the  temperature  range  of  the  flow  reactor,  and  second,  the 
high  convective  flows  of  the  flow  reactor  will  spread  the  reaction  out  over  approximately  a 
meter  in  length,  thus  providing  excellent  spatial  resolution  with  which  to  obtain  accurate 
kinetic  data.  The  model  analysis  of  RDX  flames  by  Melius  also  shows  that  the  reaction 
paths  follow  the  hierarchy  as  shown  in  FIGURE  18.  In  particular,  the  C0A320/N02  reac¬ 
tion  mechamsm  is  a  submechanism  of  both  the  CH2O/NO2  and  the  HCN7NO2  systems, 
and  the  H2/NO2  reaction  is  a  submechanism  of  the  CO/H2O/NO2  system. 

Preliminary  expenmental  results  by  Parr  (1989)  on  polycychc  mtramine  flames  has  indi¬ 
cated  that  the  interaction  between  an  energetic  binder  and  the  polycyclic  nitramine  signifi¬ 
cantly  affects  Ignition  and  combustion  kinetics,  much  more  than  with  HMX. 
Consequently,  an  understanding  of  the  gas-phase  kinetics  of  binder  intermediates,  such  as 
C2H4  and  CH4.  is  also  of  importance.  Assuming  that  the  dominant  route  for  CH3  radicals 
IS  reaction  w  ith  NO2  and  not  itself,  the  methane  reaction  may  be  considered  a  submecha¬ 
nism  of  C2H4.  However,  under  fuel-nch  conditions,  methyl  radicals  can  be  expected  to  re¬ 
combine  to  form  ethane,  and  hence  the  hierarchical  structure  of  FIGURE  18  would  have  to 
be  modified  accordingly. 

The  nitrogen  oxidizers  typically  follow  the  path  of  NO2  ->  NO  ->  N20->  N2.  Consequent¬ 
ly.  both  NO  and  N2O  as  oxidizers  are  submechanisms  of  the  systems  reponed  in  FIGURE 
18. 

A  proposed  reaction  gnd  for  RDX  combustion  is  presented  in  RGURE  20  (Alexander  et 
ai.  1989).  There  are  over  700  possible  interactions  between  the  proposed  species.  .Many 
of  these  interactions  may  be  discarded  based  on  thermochemical  and  structural  consider¬ 
ations.  Note  that  the  reaction  systems  presented  in  FIGURE  18  provide  the  necessary  build¬ 
ing  blocks  to  systematically  develop  and  validate  the  RDX  reacuon  mechanism.  With  the 
exception  of  the  last  block,  these  reaction  systems  are  also  the  building  blocks  necessary 
for  consuuicting  the  reaction  mechanisms  of  the  polycyclic  mtramines  of  interest  here. 

.■\s  indicated  in  FIGURE  18,  the  smallest  submechanism  is  the  H2/NO2  reaction  system, 
w  hich  IS  show  n  in  FIGURE  20  to  consist  of  the  well  known  H2/O2  system  and  species  and 


FIGURE  18.  Hierarchical  Development  and  Validation  of  Gas 
Phase  Reaction  Mechanisms 


reactions  important  to  NO^  formation  and  destruction.  Hence  a  sound  foundation  already 
exists  on  which  to  build  larger  mecharusms. 

Using  the  HCN/NO2  reaction  mechanism  ofThome  and  Melius  (1989)  and  retaining  only 
those  species  and  reactions  listed  in  the  H2/O2  and  NO^  formation  and  destruction  mecha¬ 
nism  of  FIGURE  20.  a  series  of  hypothetical  calculations  were  performed  on  H2/NO2  to 
simulate  expenments  in  flow'  reactors,  shock  tubes  and  premixed  flames. 

The  results  for  an  equal  molar  mixture  of  H2  and  NO2  in  nitrogen  (flow  reactor.  T,  =  1200 
K.  1  atm),  in  argon  (shock  tube.  T,  =  2200  P  =  1  atm),  and  without  diluent  (premixed 
flame.  T,  =  298  K.  P  =  0.1  atm)  are  shown  in  FIGURE  21.  The  model  predicts  that  expen- 
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FIGURE  19.  Simulated  Structure  of  RDX  Combustion  (from 
Melius,  1988). 


menis  are  possible  in  each  device,  which  will  ultimately  lead  to  validation  of  the  model. 
Compare  the  spatial  dependence  of  the  flow  reactor  and  premixed  flame.  The  results  of  a 
sensitivity  analysis  for  the  three  experiments  are  presented  in  FIGURE  22.  Note  that  the 
combination  of  reactions  important  to  the  shock  tube  and  flow  reactor  are  the  same  as 
those  imponant  to  the  flame.  However,  for  the  flame,  diffusivities  of  radicals  and  energy 
(although  not  l.sted  in  RGURE  22)  were  equally  important  parameters  in  the  model,  and 
hence  model  validation  in  the  flame  includes  both  transpon  and  chemistry.  The  cmcled  re¬ 
actions  in  FlGLT.F  22  indicate  the  weU  studied  H2/O2  reactions.  The  other  reactions  are 
from  the  NO^  fonnation  and  destruction  mechamsm.  some  of  which  have  rate  constants 
not  as  well  known. 

In  FIGURE  23.  sensitivity  analysis  results  are  presented  for  HCN/N02  mixtures,  obtained 
from  corresponding  kinetic  calculations  of  flow  reactor,  shock  tube  and  premixed  flame 
expenments.  The  circled  reactions  include  all  the  reactions  which  were  important  in  the 
H2/NO2  studies.  Note  again  that  the  reactions  imponant  to  the  flame  are  the  same  as  those 
imponant  to  the  combined  shock  tube  and  flow  reactor  studies.  Also  note  that  the  total 
number  of  important  reactions  to  both  HCN/NO2  and  H2/NO2  mixtures  is  small  compared 
to  the  total  number  of  possible  kinetic  interactions  as  described  in  FIGURE  20. 

Lastly.  FIGURE  24  presents  the  reactions  which  govern  the  RDX  bum  rate  (.Melius.  1988). 
The  circled  reactions  mclude  all  those  found  important  to  the  HCNA'02  and  H2/N'02  mix¬ 
tures.  Note  that  with  the  exception  of  N+OH  =  NO+H  (which  would  likely  show  up  as  im¬ 
portant  in  other  proposed  mixtures  to  study,  e.g..  H2/NO  mixtures),  the  only  unknown 


FIGURE  20.  Reaction  Grid  for  RDX  Decomposition 


reactions  of  importance  to  the  RDX  bum  rate  are  the  initial  decomposition  reactions  of 
RDX. 

Although  the  results  of  FIGURES  22-24  are  preliminary  and  not  comprehensive,  they  do 
emphasize  the  imponance  of  systematic  construction  of  reaction  mechanisms  and  the  im¬ 
portance  of  mechaiusm  development  via  kinetic  data  from  different  types  of  experiments. 

In  summary.  nGLT?E  18  reports  the  hierarchical  structure  of  the  proposed  research,  stan- 
ing  with  experiments  and  model  development  and  validation  for  H2  systems  and  then  sub¬ 
sequently  proceeding  towards  C2H4  and  HCN.  Several  of  these  systems,  either  with  N02 
or  \20  as  oxidizers,  have  been  studied.  However,  most  of  the  recent  studies  have  been 
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FIGURE  21,  Simulated  Experiments  with  H^'NOi  .Mixtures 


conducted  in  premixed  flame  configurations  (see  for  example.  Branch  ei  al..  1989a  and 
1989b,  Cattolica,  R.  et  ai.  1982.  Miller  and  Bowman.  1989).  Even  less  work  has  been 
conducted  on  these  mixtures  in  shock  tubes  (see  for  example  Lin,  1989).  At  present,  flow 
reactor  data  on  the  kinetics  of  these  mixtures  is  virtually  non-existent,  as  is  gas-phase  ki¬ 
netic  data  above  one  atmosphere  pressure.  Such  data  will  be  of  significant  importance  to 
the  success  of  the  overall  energetic  materials  research  program. 


N02+0-N0+02 


N02+0-N0+02 


N02+M-N0+0+M 


N02-t•M-N0■t■0-^M 

C0H-t-H2-H20t-H 

)  C0H+0H.H20-^0  ^ 

0+H2.0H+H  ) 

N02-^M=NO-^0■^M 

N02-v0-N0-f02 

OH-fNO-^M=HONO-vM 

Ch-^02=0-^0H 

)  C  H-t.02=0-^0H  ) 

H0N0-KH=H2-t-N02 

N02.kH-N0-^0H 

N02-t.N02=2N0-^02 

HNO-vOH-NO-vH20 

H02.^N0-N02-^0H 

C  0-t.H2=H»0H  ) 

H0N0-v0H=N02-t-H20 

N0-t.H02=HN0-^02 

N0-fH02=HN0->-02 

C  H02-.OH=H20*02  ) 

^H02-f0H=H20-f02  ) 

HN0-f0H=N0-^H20 

'0H+0H=H20-f0  } 

N02-^H=NO*OH 

N02-I-H-N0-I-0H 

N02-VN02-2NO-.-02 

The  sensilivities  are  evaluated  with  respect  to  the  H2  mass  traction  protiie 


FIGURE  22.  Reactions  with  High  Sensitivity  in  H2  N02  Mixtures 


3.3  Reaction  Mechanism  for  H2/N02/N2^  Kinetics 

The  initial  H2/NO2/N2O  reaction  mechanism  has  been  constructed  by  staning  with  a  com¬ 
prehensive  reaction  mechanism  for  H2/O2  kinetics.  TABLE  1  contains  the  standard  enthalp¬ 
ies  of  formation  entropies  and  specific  heats  for  the  atomic  and  molecular  species 

considtu^d  in  the  H2/O2  submecharusm.  All  of  the  data  are  from  the  JAN  AF  Tables  (1985) 
with  the  exception  of  the  heat  of  formation  for  HO2.  which  is  from  Hills  and  Howard 
(1984).  The  temperature  dependencies  of  the  data  are  stored  as  polynomial  fits  in  the  for¬ 
mat  of  the  NASA  chemical  equilibrium  program  (Gordon  and  McBride.  1971).  The  poly¬ 
nomial  coefficients  for  all  species  except  HO2  are  from  Kee  et  al.  (1987).  The  polynomial 
coefficients  for  HO2  were  obtained  using  the  THERM  code  of  Ritter  and  Bozzelli  (198'') 


HCN+O-NCO+H 


HCN+O-NCO+H 


N20+M=N2+0+M 


C  N02+M.N0+0+M 

D 

CN+02-NC0+0 

C02+N-N0+02 

C  H+02-0H+0 

D 

HCN+O-NH+CO 


HCN+0H-CN+H20 

NH2+NO-NNH+OH 
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HOCN+H-HNCO+H 

HNCO+H-NH2+CO 

CN+N20=NC0+N2 

HCN+OH-HNCO+H 

HCN+OH-HOCN+H 


C  N02f  M-NO-t-O-t-M  ) 

NCO+M-N+CO+M 
( N02-t-0.N0+02  ) 

NC0+N0-N20+C0 
( H>02-0H-t-0  ) 

N+N0-N2+0 
CN+O-CO+N 
CN+OH=NCO+H 


HCN+O-NCO+H 

NCO+M-N+CO+M 

N+N0-N2+0 


(  i,~?+M-N0+0+M 

(  N02+0-N0+02 

(  C0+0H-C02+H 

NCN+OH=HCN+NO 
HCN+0H-CN+H20 
NH*OH=HNO+H 
C  H+0H+M=.H20>M  J 


Flow  reactor  and  shock  tube  sensitivities  were  determined  with  respect  to  the  HCN  mass  traction. 
Premixed  flame  sensitivities  are  from  Thorne  and  Melius  (27th  JANNAF  Combustion  Mrotmgi 
and  were  determined  with  respea  to  the  mass  burning  rate. 


FIGURE  23.  Reactions  with  High  Sensitivity  in  HCN  Mixtures 


The  reaction  mechanism,  consisting  of  19  forward  and  reverse  elementary  reactions  is  list 
ed  in  TABLE  2.  The  forward  rate  constants  are  all  hterature  values  (see  Yetter  et  at..  1991). 
Backward  rate  constants  are  computed  from  the  forward  rate  constants  and  equilibrium 
constants.  The  temperature  range  listed  in  the  table  covers  the  domain  over  which  the  for¬ 
ward  rate  constant  parameters  were  denved.  For  many  of  the  reactions,  this  is  also  a  good 
indication  of  the  range  of  experimental  data. 


N20+M-N2+0-fM 


TABLE  1.  Hi'O;  Thermochemical  Parameters® 


SPECIES 

AH^.298 

^298 

Cp.300 

Cp,400 

Cp.soo 

Cp.600 

Cp.800 

Lp.iooo  f^p.isoo 

H 

52.10 

27.39 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

0 

59.56 

38.47 

5.23 

5.14 

5.08 

5.05 

5.02 

5.00 

4.98 

OH 

9A2 

43.88 

7.15 

7.10 

7.07 

7.06 

7.13 

7.33 

7.8 

H2 

0.00 

31.21 

6.90 

6.96 

7.00 

7.02 

7.07 

7.21 

7.73 

O2 

0.00 

49.01 

7.01 

7.22 

7.44 

7.65 

8.07 

8.35 

8.72 

H2O 

-57.80 

45.10 

8.00 

8.23 

8.44 

8.67 

9A2 

9.87 

U.26 

HO2 

3.00 

54.43 

8J6 

8.95 

9.48 

9.96 

10.75 

11.37 

12J4 

H,0, 

-32.53 

55.66 

10.41 

11.44 

12J4 

13.11 

14.29 

15.21 

16.85 

*  units  are  kcaVmol  for  AHy  29s  caVmol-K  for  S  and  Cp 


Recombinationydissociation  reaction  rate  parameters  are  reported  for  N2  as  the  bath  gas. 
The  collision  efficiencies  of  species  other  than  N2  were  accounted  for  by  including  chap¬ 
eron  efficiencies  in  the  total  concentration  [M]  as  described  in  the  bottom  of  TABLE  2 
(W'arnatz.  1985).  In  this  mechanism,  it  is  assumed  that  third  body  efficiencies,  although 
different  among  the  species,  are  the  same  for  each  recombination/dissociation  reaction. 

The  uncertainty  factor  (t'F)  is  defined  by  k  x  UF  and  =  k/UF.  The  uncertainties 
placed  on  each  rate  constant  were  derived  from  Tsang  and  Hampson  (1986)  and  Wamatz 
(1985)  and  are  intended  to  cover  the  temperature  ranges  reported. 

The  species  and  reactions  added  to  the  H2/O2  mechanism  to  account  for  N2O  or  NO2  as  an 
oxidant  of  H2  are  reported  in  TABLES  3  and  4.  respectively.  The  species  included  were 
obtained  from  a  chemical  kinetic  data  base  for  propellant  combustion  recently  generated 
by  Tsang  and  Herron  (1990).  All  of  the  thermodynamic  data  and  the  polynomial  coeffi¬ 
cients  are  from  Kee  er  al.  (1987).  The  mechanism  includes  all  the  reactions  involving  the 
species  NO.  N02.  HNO,  HONO,  and  N2O  and  their  recommended  rate  constants  as  re¬ 
viewed  by  Tsang  and  Herron  (1990).  Two  exceptions  to  this  are  the  rate  constants  for  the 
reaction  between  NO2  and  H  and  between  N2O  and  OH.  The  reaction  between  NO2  and  H 
has  been  recently  studied  by  Ko  and  Fontijn  (1991)  over  the  temperature  range  296  -  760 
K  and  we  have  adopted  their  experimental  evaluation  here.  For  the  reaction  between  N2O 
and  OH.  Tsang  and  Herron  place  an  upper  limit  on  this  reaction  and  suggest  two  product 
channels.  The  second  endothermic  channel  forms  HNO  and  NO  (A//^  295=^6.45).  The 
rate  constant  recommended  was  8.43  x  10*^  exp(-9935/RT).  As  described  below,  the  rate 
of  this  reaction  at  KXX)  K  was  found  to  be  too  large,  when  using  only  the  first  channel 
(A///J  295='25.93).  The  rate  used  by  Miller  and  Bowman  for  this  reaction  is  about  a  factor 
of  three  lower  at  this  temperature  and  is  used  here.  The  second  channel  is  included  in  the 
mechanism,  however  the  rate  constant  is  specified  for  the  backward  reaction  and  is  from 
Thome  and  Melius  (1989). 


TABLE  2.  H,  Reaction  Mechanism 


AH298 

log(Af) 

E..f 

UF 

Tra.nge 

Hi-O-,  Chain  Reactions 

1.  H+O,  =  O+OH 

16.77 

14.28 

0.00 

16.44 

2 

962-2577 

2.  0+Hi  =  H+OH 

1.86 

4.71 

2.67 

6.29 

1.5 

297-2495 

3.  OH+H,  =  H+H,0 

-15.02 

8.33 

1.51 

3.43 

1.5 

250-2581 

4.  OH+OH  =  O+HjO 

16.87 

4.09 

2.62 

-1.88 

2.5 

250-2000 

H7-O,  D'lssociation/Recombination  Reactions 

5.  H2+M  =  H+H+M 

104.1 

19.66 

-1.40 

104.38 

3.0 

600-2000 

6.  O-t-O+M  =  O2+M 

-U9.1 

15.79 

-0.50 

0.00 

1.3 

2000-10000 

7.  O+H+M  =  OH+M 

-102.3 

18.67 

-1.00 

0.00 

10.0 

8.  H+OH+M  =  H2O+M 

-U9.2 

22.35 

-2.00 

0.00 

2. 

1000-3000 

Formation  and  Consumption  of  HOi 

9.  H+0-,+M  =  HO,+M 

-49.07 

19.83 

-1.42 

0.00 

3. 

200-2000 

10.  HO1+H  =  Hi+Oi 

-55.07 

13.82 

0.00 

2.13 

2. 

298-773 

U.  HO,+H  =  OH+OH 

-36.44 

14.23 

0.00 

0.87 

2. 

298-773 

12.  H0.+0  =  0H+0i 

-53.22 

13.24 

0.00 

-0.40 

1.2 

200-400 

13.  HO2+OH  =  H2O+O2 

-70.09 

16.16 

-1.00 

0.00 

2. 

298-1400 

Formation  and  Consumption  of  H^Oi 

14.  HO7+HO,  =  H.O2+O2 

-38.51 

12.48 

0.00 

1.39 

3. 

650-800 

15.  H70,+M  =  OH+OH+M 

51.14 

17.08 

0.00 

45.50 

2. 

700-1500 

16.  h'o,+H  =  H.O+OH 

-68.02 

13.00 

0.00 

3.59 

3. 

283-800 

17.  H,0,+H  =  H7+HO, 

-17.56 

13.68 

0.00 

7.95 

5. 

250-800 

19.  H2O2+OH  =  H2O+HO2 

-32.57 

12.85 

0.00 

1.43 

2. 

98-800 

“units  are  cm'-mole-s-kcal,  k  ~  AT"  exp(-Ej  RT) 

[M]  =  fN,]+(H]+[0]+[0H]+2.5[H3]+[0J+12[H.0]+[H0J+[H>0,]+1.9[C0]-^3.8[C0,]+[HC0] 

In  addition,  a  few  reactions  have  been  added  to  the  present  system  that  were  not  reviewed 
by  Tsang  and  Herron,  w  hich  appear  in  mechanisms  developed  by  Miller  and  Bowman  for 
predicting  NO^  formation  and  destruction  and  by  Thome  and  Melius  for  predicting  HCN/ 
NO2  kinetics. 

The  pressure  dependence  of  unimolecular  reactions  is  based  on  the  relation  of  Troe 
(1977).  The  second  order  rate  constant,  k^,  the  first  order  rate  constant,  k^^f.  and  the  broad¬ 
ening  factor.  Fj..  are  also  from  Tsang  and  Herron, 

The  formation  of  NH  in  H2/N20  systems  has  been  suggested  by  Danby  and  Hinshelwood 
(1940)  and  later  by  Baldwin  etal.  (1975).  Baldwin  etal,  who  observed  NOproducuon 
near  870  K  in  bulb  experiments  attributed  the  NO  to  a  second  channel  of  the  H  +  N2O  re¬ 
action 


H  -I-  N^O  ->  SH  +  .VO  AH2q8=  34. 1  kcal/mol 

H  +  .V,0  A  ,  +  OH  AH2q8=-62.2  kcal/moie 


The  initial  mechanism  did  not  include  the  endothermic  channel  because  NO  was  not  ob¬ 
served  in  the  H2/N2O  experiments.  How  ever,  to  check  the  importance  of  this  channel  and 
the  subsequent  reactions  which  consume  NH.  reactions  53  -  73  in  TABLE  4  were  includ¬ 
ed.  which  are  from  Miller  and  Bowman  and  Thorne  and  Melius.  Reactions  involving  HN2. 
HN3,  and  NNH  were  not  included  because  there  is  little  likelihood  that  these  species  play 
any  significant  role  in  the  H2/N2O  and  H2/NO2  systems  of  this  work. 


TABLE  3.  Nitrogen  Species  and  Thermochemical  Parameters*  added  to  the  H2/O2  System 


SPECIES 

4Hy^,29S 

S»8 

Cpjoo 

Cp.400 

Cp.500 

Cp.600 

Cp.800 

Cp.lOOO  f^p.1500 

NO 

21.58 

50.35 

7.U 

7.19 

7J1 

7.45 

7.82 

8.14 

8.54 

NO, 

7.91 

57.34 

8.83 

9.64 

10J3 

10.93 

U.89 

12.49 

13.17 

NjO 

19.61 

52.55 

9.27 

10.18 

10.94 

11.56 

12jl 

13.12 

13.94 

HNO 

23.80 

52.73 

8.26 

8.84 

9J6 

9.84 

10.76 

11.48 

12.49 

MONO 

-18J4 

59.59 

10.88 

12.27 

13_39 

14.31 

15.67 

16.56 

17.89 

MONO, 

-30.18 

63.17 

12.34 

14.70 

16.61) 

18.12 

20.30 

21.70 

23.45 

NO., 

17.00 

60.37 

11.31 

13.32 

14.90 

16.10 

ns9 

18.27 

19.10 

N 

112.96 

36.61 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

NH 

85.21 

43.30 

6.96 

6.98 

7.00 

7.04 

7.21 

7.48 

8.06 

*  units  are  kcal/mol  for  AH 

f.29S  CQL" 

mol-K  for  S  and  Cp 

T.\BLE  4.  NO;  N;0  Reaction  Mechanism  Added  to  the  H;  O;  System 


Af 

"f 

Ea.f 

Ref. 

NO  Reactions 

20.  N0+M=N+0+M 

0.14.'E4.16 

0.000 

148400.000 

Tsang  and  Herron  (1990) 

21.  N0+0=N02 

kiin/)  0.130E+16 

-0.750 

0.000 

Tsang  and  Herron  (1990) 

k(0,N,)  0.472E+25 

-2.870 

1552.000 

FciNN)  =0.95-10  *r 

22.  N0+0=02+N 

0.181E+10 

1.000 

38750.000 

Tsang  and  Herron  (1990) 

23.  NO+H=HNO 

kdnf)  0.152E+16 

-0.410 

0.000 

Tsang  and  Herron  (1990) 

k(0,N2)  0.896E+20 

-1.320 

735.200 

F,(N2)=0.82 

24.  NO+H=N+OH 

0.169E+15 

0.000 

48800.000 

Tsang  and  Herron  (1990) 

25.  NO+OH=HONO 

k(inf)  0.199E+13 

-0.050 

-7.213 

Tsang  and  Herron  (1990) 

klO.Nj)  0.508E+24 

-2.510 

-67.560 

Fr(N2)  =  0.62 

26.  2N0rN2+02 

0.130E+15 

0.000 

75630.000 

Thorne  and  Melius  (1989) 

27.  N0+H02=HN0+02 

0.200E*12 

0.000 

2000.000 

Thorne  and  Melius  (1989) 

NO,  Reactions 

28.  N02+H2=H0N0+H 

0.24IE+14 

0.000 

28810.000 

Tsang  and  Herron  (1990) 

29.  N02*0=02*N0 

0.^91  E*13 

0.000 

-2.^8.400 

Tsang  and  Herron  '  1990) 

30.  N02+0=N0, 

k(inf)  0.133E+14 

0.000 

0.000 

Tsang  and  Herron  (1990) 

k(0,N2)  0.149E+29 

-4.080 

2468.000 

Ft(N2)  =0.79-1.8xl0-*T 

31.  N02+H=N0+0H 

0.129E+15 

0.000 

361.600 

Ko  and  Fontijn  (1991) 

32.  N02+0H=H0N02 

k(in/)  0.241E+14 

0.000 

0.000 

Tsang  and  Herron  (1990) 

k(0,N2)  0.642E+33 

-5.490 

2351.000 

F^CN'i)  =  0.725-2.5x10^ 

33.  N02+0H=H02+N0 

0.181E+14 

0.000 

6676.000 

Tsang  and  Herron  (1990) 

34.  N02+N0=N20+02 

Q.mE*13 

0.000 

60000.000 

Thorne  and  Melius  (1989) 

35.  2N02=N03+N0 

0.964E4^10 

0.730 

20920.000 

Tsang  and  Herron  (1990) 

36.  2N02=2N0+02 

0.163E+13 

0.000 

26120.000 

Tsang  and  Herron  (1990) 

NjO  Reactions 

37.  N20=N2+0 

k(in/|  0.130E+12 

-0.000 

59610.000 

Tsang  and  Herron  (1990) 

k(0,N,)  0.723E+18 

-0.730 

62780.000 

F^(N2)  =  1.167-1.25x10  ^ 

38.  N204.0=N2+02 

0.102E+15 

0.000 

28020  000 

Tsang  and  Herron  (199G) 

39.  N20+0=2N0 

0.663E+14 

0.000 

26630.000 

Tsang  and  Herron  (1990) 

40.  N20+H=N2+0H 

0.964E+14 

0.000 

15100.000 

Tsang  and  Herron  (1990) 

41.  N20+0H=H02+N2 

0.200E+13 

0.000 

10000.000 

Miller  &  Bowman  (1989) 

42.  N20+N0=N2+N02 

O.lOOE+15 

0.000 

49675.000 

Thorne  and  Melius  (1989) 

HNO  Reactions 

43.  HN0+0=0H*N0 

0.361E+14 

0.000 

0.000 

Tsang  and  Herron  (1990) 

44.  HNO+0=N02+H 

0.500E+11 

0.500 

2000.000 

Thorne  and  Melius  (1989) 

4.'.  HN0+H=H2+N0 

0.181E+14 

0.000 

993.500 

Tsang  and  Herron  (1990) 

46.  H\0.^0H=H20+N0 

0.482E+14 

0.000 

993.500 

Tsang  and  Herron  (1990) 

47.  HN0+N0=N20+0H 

0.200E+13 

0.000 

26000.000 

Thorne  and  Melius  (1989) 

48.  HN0+N02=H0N0+N0  0.602E+12 

0.000 

1987.000 

Tsang  and  Herron  (1990) 

49.  2HN0=H20+N20 

0.843E+09 

0.000 

3102.000 

Tsang  and  Herron  (1990) 

MONO  Reactions 

50.  H0N0+0=0+N02 

0.120E+14 

0.000 

5961.000 

Tsang  and  Herron  (1990) 

51.  H0N0+0H=H20+N02  0.126E+n 

1.000 

135.100 

Tsang  and  Herron  (1990) 

52.  H0N0+HN0=H20+2N0  O.lOOE+13 

0.000 

40000.000 

Miller  &  Bowman  (1989) 

NH  Reactions 

S3.  NH+02=HN0*0 

O.lOOE+14 

0.000 

12000.000 

Miller  &  Bowman  (1989) 

54.  NH+0,=N0+0H 

0.760E*U 

0.000 

1530.000 

Miller  &  Bowman  (1989) 

55.  N+HO,=NH+Oj 

O.lOOE+14 

0.000 

2000.000 

Miller  &  Bowman  (1989) 

56.  NH+0=N0+H 

0.200E+I4 

0.000 

0.000 

Miller  &  Bowman  (1989) 

57.  NH+HrN+H, 

0.300E+14 

0.000 

0.000 

Miller  &  Bowman  (1989) 

58.  N+0H=NH+0 

0.380E+14 

0.000 

0.000 

Miller  &  Bowman  (1989) 

59.  .NH+OH=HNO.^H 

0.200E*14 

0.000 

0.000 

Miller  &  Bowman  (1989) 

60.  NH+0H=N*H,0 

0.500E+12 

0.500 

2000.000 

Miller  &  Bowman  (1989) 

61.  NH+N=N,+H 

O.lOOE+15 

0.000 

0.000 

Miller  &  Bowman  (1989) 

62.  NH+NOrNiO+H 

0.240E+16 

-0.800 

0.000 

Miller  &  Bowman  (1989) 

63.  N*HNO=NH+NO 

O.lOOE+14 

0.000 

2000.000 

Miller  &  Bowman  (1989) 

64.  NH+NO.rNO+HNO 

O.lOOE+12 

0.500 

4000.000 

Miller  &  Bowman  (1989) 

% 

% 


% 


N  Reactions 


65.  N+N0=N2-‘-0 

0.327E+13 

0.300 

0.000 

Miller  & 

Bowman  (1989) 

66.  N+N20=N,-^N0 

O.lOOE+14 

0.000 

19870.000 

.Miller  & 

Bowman  (1989) 

67.  N+N02=2N0 

0.400E+13 

0.000 

0.000 

Miller  & 

Bowman  (1989) 

68.  N+N02=N20+0 

0.500E+13 

0.000 

0.000 

Miller  & 

BowTTian  (1989) 

69.  N+N02=N2+02 

O.lOOE+13 

0.000 

0.000 

Miller  & 

Bowrman  (1989) 

70.  N+HN0=N20+H 

0.500E+U 

0.500 

3000.000 

Miller  & 

Bowman  (1989) 

71.  N+H02=N0+0H 

O.IOOE.^14 

0.000 

2000.000 

Miller  & 

Bowman  (1989) 

•units  are  cm^-moU-s-cal,  k  -  AT"  exp(-E,/RT) 


4.0  N2O  DECOMPOSITION  AND  REACTION  WITH  H2 

4.1  Introduction 

A  number  of  studies  exist  for  the  reaction  between  H2  and  N2O.  However,  most  of  the  re¬ 
cent  studies  are  for  high  temperatures,  typically  from  expenments  conducted  in  flames  or 
shock  tubes.  Flames  studies  include  the  low  pressure  (40  lonr)  lean  H2/N2O  flames  of  Bal- 
akhnine  et  al.  (1977).  the  atmospheric  pressure  stoichiometric  H2/N2O  flame  of  Cattolica 
et  al.  (1982)  and  the  flame  study  of  Fine  and  Evans  (1963).  Recent  shock  tube  studies  of 
H2/N2O  mixtures  have  often  been  performed  to  measure  the  rates  of  isolated  reactions. 
For  example.  Hidaka  et  al  ( 1985)  studied  the  rate  of  N20  +  H  and  Pamidimukkaln  and 
Skinner  (1982)  studied  the  rate  of  H2+O.  Henrici  and  Bauer  (1964)  used  a  shock  tube  to 
study  this  system  at  1700-2600  K  and  obtained  a  value  for  the  rate  of  K-t-N^O 

In  a  temperature  range  slightly  below  that  of  this  study.  Baldwin  et  al.  (19“5)  have  studied 
the  H2/N2O  reaction  in  a  bulb  reactor  over  the  temperature  range  of  810  -  870  K.  Their 
work  suggested  the  imponance  of  NH  and  HNO  in  this  system.  These  reactions  were  con¬ 
sidered  to  lead  to  termination  in  the  reaction.  Further,  wall  reacuons  could  not  be  neglect¬ 
ed. 

The  present  studies  have  been  conducted  over  a  temperature  range  too  high  for  bulb  ex¬ 
periments  and  too  low  for  shock  tubes  (950  -  1073  K)  in  a  flow  reactor  where  wall  reac¬ 
tions  are  minimized.  Detailed  kinetic  modeling  and  sensitivity  analyses  are  also  performed 
to  thoroughly  investigate  the  model. 

4.2  Results 

TABLE  5  reports  the  initial  experimental  conditions  of  the  systems  studied  for  N2O  de¬ 
composition  and  Its  reaction  with  H2.  An  example  of  the  species  profiles  for  N2O  decom¬ 
position  IS  reported  in  FIGURE  25  for  the  condiuons  of  EXP  2.  Approximately  14%  of  the 
N^O  decomposes  during  an  experimental  reaction  time  of  0.75  s.  During  this  penod.  ap¬ 
proximately  490 ppm  of  O2.  50  ppm  of  NO.  and  4 ppm  of  NO2  are  formed.  For  the  same 
reaction  time,  the  results  of  experiment  EXP  1  showed  approximately  6%  N20  consump¬ 
tion  and  formation  of  220 ppm  of  O2.  20 ppm  of  NO  and  2  ppm  of  NO2.  Hence,  the  rate  of 
reaction  is  nearls  proportional  to  pressure. 


TABLE  5.  Experimental  Conditions  of  S2O  Experiments 


EXP 

TX 

P/atm 

1 

1075 

32 

2 

1096 

6.0 

3 

950 

3.0 

4 

1014 

3.0 

5 

1020 

3.0 

6 

1073 

1.0 

7 

1073 

1.0 

8 

1025 

3.0 

Xn20 

Xh2 

<b 

0.008 

0.0 

0.078 

0.0 

0.009 

0.0097 

1.08 

0.01 

0.0072 

0.72 

0.0074 

0.0052 

0.70 

0.0074 

0.0074 

1.0 

0.0065 

0.0065 

1.0 

0.0070 

0.070 

1.0 

The  results  suggest  that  the  rate  limiting  step,  as  expected,  is  the  thermal  decomposition  of 
N2O  via  N2O  +  M  ->  N2  +  O  +  M  and  that  the  branching  ratio  of  N2O  +  O  yielding  N2  + 
O2  vs.  NO  +  NO  is  20  to  1 .  However,  according  to  TABLE  4.  the  ratio  of  the  rate  constants 
for  reaction  38  to  that  of  reaction  39  at  1096  K  is  0.81.  This  mechamsm  implies  NO  is  be¬ 
ing  consumed  or  O2  is  being  formed  from  other  reactions  than  N2O  -t-  O.  or  that  the 
branching  ratio  given  in  TABLE  4  is  incorrect. 

The  corresponding  modeling  calculations  are  shown  in  FIGURE  26.  The  model  does  not 
predict  the  fast  reaction  for  times  less  than  0.1  s.  but  generally  agrees  with  the  decomposi¬ 
tion  rate  for  times  greater  than  0.1 5.  The  model  predicts  about  8.5%  consumption  of  N2O. 
formation  of  230  ppm  of  O2.  84  ppm  of  NO  and  5  ppm  of  NO2  for  a  reaction  time  of  0.75 
5. 

The  corresponding  sensitivity  analysis  results  for  the  gradients  dY\-2o/^^"^net 

sented  in  FIGURE  27.  The  most  sensitive  reactions  in  decreasing  order  of  imponance  are. 


A'jO  +  M  N2  +  O  +  M  37 

2  G  +  O  ^  A  T  +  O  T  38 

NO  +  0  -k-M  +  M  21 

A'jO  ^-Oi^NO^NO  39 

NO2  ■^0^r^02  +  ^'0  29 


The  numbers  above  and  in  the  figures  denote  the  reaction  numbers  as  listed  in  TABLES  2 
and  4.  The  important  parameters  of  the  reactions  above  are  the  forward  rate  constants. 

The  effect  of  O2  and  H2O  impurities  in  the  N2  carrier  gas  are  illustrated  in  FIGURE  28  for 
concentrauons  of  250  ppm  each  added  separately  to  the  N2O  system  of  RGURE  26.  The 
mechamsm  used  here  only  included  reactions  1  -  52  of  TABLES  2  and  3.  The  model  results 
indicate  that  O2  impurities  do  not  have  a  significant  affect  on  the  decomposition  rate.  In 
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Concentration  (ppm) 


N20  Decomposition 

6atm,  1096K 


FIGU  RE  25.  Experimental  concentration  profiles  of  N^O 
decomposition.  Initial  conditions  are  reported  in  TABLE  5  for  EXP 
2 


contrast,  water  vapor  has  an  appreciable  affect  on  the  rate  of  N2O  decomposition.  Sensi¬ 
tivity  analysis  results  indicate  that  in  addition  to  reactions  33  and  39  above,  the  reactions 


HO2  +  A’O  AOj  +  OH  33 

HO2  +  O  i-^OH  +  O2  12 

HO2  +  OH  H^O  +  O2  13 

0H  +  .\20<r^H02  +  S2  41 


.01  — 


0  .2  .4  .6  .8  1 


0  .2  4  .6  ,6  1 

l/s 


FIGURE  26.  Calculated  mass  fraction  profiles  of  NiO 
decomposition.  Initial  conditions  arc-  reported  in  TABLE  5  for  EXP 
2 


H  +  0^<r4  0H  +  0  ] 

O  +  H^O^OH  +  OH  4 

are  imponant  when  H2O  is  present.  The  addition  of  the  NH  reactions  and  N  reactions  had 
no  affect  on  the  water  catalyzed  N2O  decomposition  reaction  because  H  atoms  are  not 
formed  in  significant  quantities  in  this  system. 

An  example  of  the  kinetics  of  N2O  reacting  with  H2  is  given  in  FIGURE  29  for  the  condi¬ 
tions  of  EXP  6.  Dunng  the  reaction  period  of  0. 18  s,  50%  of  the  N2O  and  H2  are  con¬ 
sumed.  No  NO  was  observed  in  these  expenments.  The  corresponding  calculated 
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FIGURE  27.  Sensitivity  gradients  dy\-,Q'dlnk^  for  the  response 
of  the  N^O  mass  fraction  to  relative  v'anations  in  trie  net  rate 
constants.  The  initial  conditions  are  the  same  as  those  of  FIGURE 
26. 


temperature  and  mass  fraction  profiles  for  N2O.  H2O.  H2.  H,  O.  OH.  HO2,  NO.  and  O2  for 
the  stoichiometric  H2/N20  system  reacting  at  1073  K  and  1  atm  are  presented  in  FIGURE 
30.  Note  the  formation  of  molecular  oxygen  and  its  consumption  afia  approximately  60% 
of  the  imtial  H2  and  N2O  have  been  consumed.  The  model  predicts  50%  consumption  of 
N2O  in  about  the  same  period  of  time  as  the  expenment. 

The  sensitivity  gradients.  stoichiometnc  H2/ 

N2O  system  reacting  at  1073  K  and  1  atn\  are  presented  in  FIGURE  31.  The  correspond¬ 
ing  gradients  for  both  species  have  opposite  signs,  but  are  of  similar  shape.  Both  species 
are  sensitive  to  the  same  reactions.  In  decreasing  order  of  importance,  these  reactions  are; 


FIGL'RE  28.  (a)  and  (b):  .Mass  fraction  profiles  for  NiO 
decomposition  in  the  present  of  250  ppm  Oi.  (c)  and  (d):  Mass 
fraction  profiles  for  NiO  decomposition  in  the  present  of  250  ppm 
HiO.  The  initial  conditions  are  the  same  as  those  of  FIGURE  26, 
but  with  the  addition  of  250  ppm  or  250  ppm  H;0. 


N2O  M  Nj-^r  O  +  M 
+  H^O  +  H 
N^O  +  OH  + 

H  +  H02<r^  H2  +  O2 

H  +  H02<r^OH  +  OH 


N.O  +H  ir^N.  +  OH 


Hydrogen  /  Nitrous  Oxide 

Phi=1, 1  Atm,  1073K,  (VPFR) 


*-  N20  (CO  meter)  H20  (FTIR)  0-bal  (CO  meter) 


FIGURE  29.  Experimental  concentration  profiles  for  the  reaction  of 
Hi  and  NiO  with  the  conditions  of  EXP  6  of  TABLE  5. 


H  +  Oj^^OH+O  1 

Using  the  upper  limit  for  the  reaction  between  OH  and  N2O  recommended  by  Tsang  and 
Hampson  and  assuming  only  the  product  channel  to  form  HO2  and  N2.  the  overall  reaction 
IS  approximately  a  factor  of  two  slower  than  the  expenment.  Note  that  among  the  impor¬ 
tant  reacoons,  this  reaction  is  the  one  with  the  largest  uncertainty. 

In  FIGURE  32,  the  gradients  of  BYj^,2o/dlnk^  and,  kinetic  sys¬ 

tem  as  above  are  reponed.  These  results  illustrate  that  for  all  the  reactions  listed  above,  the 
fopx  ard  rate  constant  is  the  imponant  parameter;  only  for  reaction  3  is  the  backward  rate 


FIGURE  30.  Calculated  temperature  and  mass  fraction  profiles  for 
the  reaction  of  H  7  and  NjO  with  the  conditions  of  EXP  6  of  T.\BLE 
5. 


constant  somewhat  important.  For  this  system,  none  of  the  equilibrium  constants  play  a 
sigmficant  role. 

The  importance  of  oxygen  in  this  system  is  illustrated  in  FIGURE  33.  In  this  figure,  the 
N2O  profiles  from  stoichiometric  mixtures  of  H2/N2O  reacting  in  N2  at  1073  K  and  1  atm 
are  presented  for  the  addition  of  trace  amounts  of  O2  ranging  from  0  ppm  to  500  ppm.  Af¬ 
ter  the  addition  of  about  50  ppm,  the  reaction  character  changes  complnely  to  a  branched 
chain  system, 

The  sensitivity  gradients  for  the  system  with  250 ppm  O2  is  reported  in 

FIGURE  34.  In  decreasing  order  of  imponance.  the  most  sensitive  reactions  are  1 1 ).  (3). 


(41),  (10),  ( 1 1).  (37),  (9),  and  (39).  Although  the  order  of  most  sensitive  reactions  has 
changed  relative  to  the  system  without  molecular  oxygen,  the  list  contains  the  same  set  of 
reactions. 

The  importance  of  NH  in  this  mechanism  is  illustrated  in  FIGURE  35.  The  calculated 
overall  rate  of  reaction  is  observed  to  decrease  by  approximately  25%.  The  effect  is  con¬ 
sistent  with  that  found  by  Baldwin  etal.  A  concentration  of  approximately  \0ppm  NO  is 
predicted  at  the  largest  extent  of  reaction  studied  experimentally  (whereas  for  the  calcula¬ 
tions  without  NH,  the  concentration  of  NO  predicted  was  less  than  0. 1  ppm).  Although 
\0ppm  NO  is  clearly  within  the  detectability  limits  of  the  FTIR,  NO  was  not  observed  in 
I  the  experiments. 
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nGL'RE  32.  Sensitivity  gradients  d)\-^idlnk.  and  d}’. 
for  the  conditions  of  FIGURE 


Sensitivity  analysis  results  reveal  that  the  inclusion  of  NH  to  the  mechanism  add  four 
more  reactions  to  those  found  important  above. 

A';C>  +  Hir^NH  +  SO  62 

H  +  SO  +  M<r^  HSO  +  M  23 

SO  +  H  <r^S  +  OH  24 

+  A  6?  A  2  +  O 
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FIGURE  33.  Calculated  NSO  mass  fraction  profiles  for  the  initial 
conditions  of  EXP  6  of  TA6LE  5  with  various  quantities  of  O, 
added.  Oj  concentrations  include  0,  10,  20,  30,  40,  50, 100,  and  250 
PP>»  _ 


Finally,  reactions  involving  NH2.  NH3,  and  NNH  in  the  NO^  mechanism  of  Melius  and 
Bowman  were  added  to  this  system  and  found  to  have  no  effect. 

4.3  Discussion  and  Conclusions 

The  results  of  N2O  decomposition  show  that  the  model  and  experiment  are  in  general 
agreement.  However,  further  experiments  are  needed  to  validate  the  model.  The  experi¬ 
ments  need  to  be  extended  to  longer  extents  of  reaction  and  higher  pressures.  Trace  quan- 
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FIGURE  34.  Sensitivity  gradients  dY .^-.^pldlnk for  the 
conditions  of  E.XP  6  with  250  ppm  Oj 

titles  of  water  vapor  have  been  shown  numencally  to  significantly  affect  the  rate  of  N2O 

decomposition  due  to  the  importance  of  OH  radicals  reacting  with  N20.  The  OH  radicals 

are  formed  from  the  branching  reaction  between  H2O  and  0-atoms,  the  latter  which  are 

formed  from  the  thermal  decomposition  reaction.  The  products  of  the  reaction  between 

OH  and  N2O  are  N2  and  HO2.  The  HO2  continues  the  chain  by  reacting  with  the  NO  ^ 

formed  from  the  reaction  of  0-atoms  with  N2O.  The  products  of  the  HO2  +  NO  reaction 

are  NO2  and  OH.  in  which  OH  is  reformed  allowing  the  chain  to  start  over  again.  Note 

that  H-atoms  are  not  formed  m  significant  quantities  in  this  system  and  hence  the  reactions 

of  H-i-N^O  are  not  important.  Further,  without  H-atoms.  the  molecular  oxygen  which  is 

formed  from  N2O  and  NO2  reacting  with  O-atoms  does  not  lead  to  branching,  as  it  can  in  a 

systems  with  H2  present  instead  of  H2O.  Experiments  conducted  with  various  amounts  of  ^ 
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FIGURE  35.  Calculated  temperature  and  mass  fraction  profiles  for 
the  reaction  of  Hi  and  N-.0  for  the  conditions  of  EXP  6  of  TABLE  5 
with  reactions  51-73  of  TABLE  4  added  to  the  mechanism. 


water  vapor  are  therefore  important  and  will  be  conducted.  Trace  amounts  of  O2  were 
found  numerically  to  have  little  affect  the  decomposition  rate. 

The  results  of  the  H2/N2O  reaction  again  show  the  model  to  generally  predict  the  experi¬ 
ments.  The  results  presented  here  suggest  that  at  the  present  conditions,  the  reaction  of  H 
+  N2O  ->  NH  +  NO  may  not  be  imponant.  The  rate  constants  presently  available  in  the  lit¬ 
erature  suggest  that  the  inclusion  of  NH  in  the  mechanism  would  slow  the  overall  reaction 
by  approximately  259c.  However,  the  presence  of  NH  in  the  mechanism  predicts  higher 
NO  levels  than  were  measured  in  the  experiments.  For  the  H2/N2C)  systems,  the  reaction 
between  OH  and  N2O  was  again  found  to  be  important  at  our  conditions.  The  importance 
of  this  reaction  and  of  HO2  radicals  in  the  H2/N2^  system  appears  to  be  greater  than  pre- 


viousiy  observed  in  flame  studies  or  in  the  lower  temperature  study  of  Baldwin  er  ai. 
(who  did  not  consider  it  at  all).  The  results  also  show  (which  are  consistent  with  other 
studies,  e.g.,  Balakhnine  et  al.)  that  the  formation  of  molecular  oxygen  is  important  in 
pure  H2/N2O  systems  because  it  leads  to  branching  through  reaction  of  H-atoms  with  mo¬ 
lecular  oxygen.  Studies  with  addition  of  molecular  oxygen  are  therefore  important  and  can 
be  used  to  adjust  residence  times  in  the  experiments  for  various  reaction  temperatures. 

5.0  NO2  DECOMPOSITION  AND  REACTION  WITH  H2 

5.1  Introduction 

At  intermediate  temperatures,  nitrogen  dioxide  is  reduced  only  panially  by  hydrogen,  ac¬ 
cording  to  the  stoichiometry. 

H2  +  NO2-*  +  Nd 

At  higher  temperatures,  such  as  in  flames,  complete  reduction  of  nitrogen  dioxide  occurs. 

The  work  of  Ashmore  and  Levitt  (1956.  1957.  1958)  and  of  Rosser  and  Wise  (1957)  on 
the  thermal  reaction  between  hydrogen  and  nitrogen  dioxide  have  convincingly  shown 


that  the  mechanism  is  straight  chain  consisting  mainly  of. 

S’02  +  H2^H0S0-^H  28 

H  +  S02<r^OH  +  SO  31 

OH  +  H2<r^H20  +  H  3 

0H  +  N02  +  M^HS0^  +  M  32 

OH  -i-  SO  +  M  HOSO  +  M  25 


The  rate  of  reaction  in  the  intermediate  temperature  range  was  experimentally  found  to  be 
empirically  expressible  as 

dlS02]  /([AO,] 

^dt  '  [NdV-r[NO^~  ■ 

According  to  this  relationship,  the  reaction  rate  is  self-inhibiting  because  of  NO  forma¬ 
tion.  which  leads  to  chain  termination.  Steady  state  analysis  (Rosser  and  Wise,  1957)  of 
[H]  and  [OH]  of  the  above  set  of  elementary  reactions  yields 

d[N02]  _  ^28^3  [AOj]  [f/j]  ^ 

dt  [AO]  +  kj2  [AOj] 

The  lower  [H2]  reaction  order  observed  expenmentally  was  attributed  to  additional  H- 
atom  loss,  e.g.,  via  wall  termination. 


Sawyer  and  Classman  (1968)  studied  the  H2/NO2  reaction  in  a  flow  reactor  similar  to  the 
one  used  in  the  present  study  (where  wall  termination  should  be  negligible).  \ei  Sawyer 
and  Classman  found  a  H2  dependence  similar  to  that  of  Rosser  and  Wise  and  of  Ashmore 
and  Levin,  where  wall  termination  was  present.  They  suggested  additional  termination 
through  the  reaction  H  +  NO  +  M  ->  HNO  +  M  and  showed  that  by  considering  this  reac¬ 
tion  as  the  main  terminating  step  that  the  theoreucal  H2  dependence  is  lowered  to  unity. 
Thus  any  termination  reaction  removing  hydrogen  atoms  acts  to  reduce  the  order  with  re¬ 
spect  to  hydrogen. 

Anderson  et  al.  (1989)  have  recently  performed  a  detailed  modelling  study  of  the  homoge¬ 
neous  igmtion  kinetics  of  H2/NO2  mixtures.  The  two  stage  process  of  the  nitrogen  chemis¬ 
try  was  shown  to  occur  through  conversion  of  NO2  to  NO  by  reaction  with  H  atoms  in  the 
first  stage  and  conversion  of  the  NO  to  N2  by  N  atoms  or  via  the  formation  and  consump¬ 
tion  of  HNO.  NH.  and  N2O  in  the  second  stage. 

Volpom  and  Branch  (1990)  have  recently  studied  the  H2/NO2  system  in  laminar  premixed 
flames.  Agreement  between  experiment  and  numerical  model  predictions  were  generally 
in  good  agreement.  The  reaction  was  found  to  be  most  sensitive  to  the  reactions 


H  +  S02ir^0H  +  S0  31 

+  +  H  3 

O  +  H^^OH  +  H  2 

H  +  O.^OH  1 

OH  +  OHir-^O  +  H^O  4 

A’O;  +  M  4-4  ,\’0  +  0  +  M  21b 

NO2  +  0  4^  \0  +  O2  29 


In  contrast  to  intermediate  and  lower  temperatures  where  molecular  oxygen  and  O-atorn 
reactions  were  found  to  be  ummportant.  the  flame  results  show  that  these  reactions  are  of 
sigmficant  importance.  This  change  in  mechanism  results  from  the  increased  rate  of  NO2 
decomposition  at  higher  temperatures  yielding  0-atoms  which  react  with  NO2  to  produce 
02.  Because  of  the  branching  introduced  into  the  system  via  formation  of  O2.  the  reaction 
of  H  NO2  ->  NO  +  OH  acts  to  inhibit  the  overall  rate. 

5.2  Results 

TABLE  6  reports  the  initial  experimental  conditions  of  the  systems  studied  for  NO2  de¬ 
composition  and  Its  reaction  with  H2.  An  example  of  the  experimental  results  of  H2  react¬ 
ing  with  NO2  IS  reported  in  FIGURE  36  for  the  initial  conditions  of  EXP  4  listed  in 
TABLE  6. 


TABLE  6.  Experimental  Conditions  of  NOi  Experiments 


EXP 

T/K 

P/atm 

Xa: 

<I> 

1 

1073 

3 

0.0002 

0.0 

2 

953 

3 

0.0002 

0.0 

3 

853 

3 

0.0002 

0.004 

10 

4 

873 

3 

0.0005 

0.004 

4 

The  experiments  show  the  consumption  of  NO2  and  the  formation  of  H2O.  A  small  tem¬ 
perature  rise  is  also  observed.  The  formation  of  a  currently  unidentified  species  (speculat¬ 
ed  as  HNO3  or  HONO)  is  observed  to  occur  in  appreciable  amounts.  The  FTIR  spectra 
and  calibration  of  these  species  is  currently  being  investigated.  The  species  shown  in  the 
figure  IS  assumed  to  be  HNO3  and  is  calibrated  by  0-atom  balance  closure.  Further  studies 
need  to  be  conducted  to  verify  that  this  species  is  a  gas  phase  product  and  not  a  result  of 
the  sampling  process. 

Model  predictions  for  this  experiment  are  presented  in  FIGURE  37.  As  can  be  seen  from 
the  figure,  the  experiment  is  about  a  factor  of  2-3  slower  than  the  model  predicts.  Further, 
formation  of  either  HONO  or  HNO3  in  quantities  no  greater  than  4  ppm  are  predicted.  Of 
the  three  radicals.  H.  O.  and  OH.  only  OH  is  present  in  appreciable  concentrations. 

The  corresponding  sensitivity  analysis  results  for  the  gradient.  present¬ 

ed  in  FIGURE  38.  The  most  sensitive  reactions  in  decreasing  order  of  imponance  are 


0// -t- WjO  + //  3 

A  O,  +  W,  HOSO  +  H  28 

HONO  +  OH  ir^H^O  +  NO  51 

H  +  .\02‘^0H  +  .\0  31 

0H-^-^0  +  M<r^H0S0  +  M  25 

SO  +  H  +  Mir^HSO  +  M  23 


Reactions  3.  28.  31 .  and  25  are  in  agreement  with  previous  studies.  However,  the  present 
model  also  indicates  the  importance  of  reaction  51.  which  consumes  HONO.  and  of  reac¬ 
tion  23.  which  removes  H-atoms  from  the  system.  Hence,  both  of  these  reactions  inhibit 
the  reaction  rate. 

Of  the  above  reactions  with  the  greatest  sensitivity,  reactions  28  and  51  have  the  largest 
uncertainties.  The  recommended  value  for  reaction  28  is  from  Slack  and  Grillo  (1978) 
who  derived  the  rate  from  an  analysis  of  the  complex  kinetics  of  H2/O2/NO2  igmtion.  An 
uncenainty  of  a  factor  of  5  is  assigned  by  Tsang  and  Herron.  The  experiments  of  Slack  and 
Grillo  were  conducted  with  non-dilute  mixtures  in  a  shock  tube.  Significant  pressure  (fac¬ 
tor  of  3)  and  temperature  (approx.  200  K)  rises  are  predicted  from  numerical  simulations 
of  their  expenments.  Hence,  their  results  are  highly  coupled  to  the  temperature  and  pres- 


Concentration  (ppm) 


N02  +  H2  reaction 


FIGURE  36.  Experimental  concentration  profiles  for  the  reaction  of 
Hi  and  NO2  with  the  conditions  of  EXP  4  of  TABLE  6. 


sure  histones  of  the  experiments,  which  were  not  momtored.  FIGURE  39  reports  the  nu- 
mencal  results  for  the  experimental  condition  EXP  4  with  a  factor  of  5  decrease  in  the  rate 
of  reacDon  28.  The  variation  is  the  uncertainty  limit  and  produces  good  agreement  be¬ 
tween  the  expenmental  and  model  predicted  rates  of  NO2  disappearance. 

The  value  for  reaction  51  is  also  that  recommended  by  Tsang  and  Herron  and  is  derived 
from  low  temperature  data.  The  only  high  temperature  data  (1000  K)  is  from  Fifer  (1976), 
which  IS  about  a  factor  of  10  lower  than  the  value  used  here.  Tsang  and  Herron  report  an 
uncertainty  of  a  factor  of  2  on  their  low  temperature  extrapolation.  A  decrease  in  reaction 
5 1  would  not  improve  the  agreement  in  the  rates  of  NO2  disappearance  between  experi- 
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FIGURE  37.  Calculated  temperature  and  mass  fraction  profiles  for 
the  reaction  of  Hi  and  NO,  with  the  conditions  of  EXP  4  of  TABLE 
6 


meni  and  model  prediction.  However,  this  change  may  yield  more  MONO  as  a  product 
species. 

The  value  used  for  reaction  31  is  from  Ko  and  Fontijn  w  ho  note  that  this  reaction  has  been 
reponed  to  proceed  through  a  MONO*  intermediate. 
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and  suggest  for  propellant  combustion  that  the  stabilization  channel  leading  to  HONO  be 
included  in  kinetic  mechanisms.  They  estimate  that  the  stabilization  rate  is  approximately 
equal  to  the  dissociation  rate  for  [M]  approximately  10^'  cm'^.  The  stabilization  channel 
would  slow  the  reaction  rate  in  the  model  predictions  here.  However,  the  pressures  of  the 
present  experiments  produce  values  of  [M]  two  orders  of  magnitude  lower  than  this  esti¬ 
mate. 

In  the  previous  section,  oxygen  addition  to  the  H2/N2O  system  was  observed  to  accelerate 
the  rate  of  reaction.  FIGURE  40  shows  the  effect  of  O2  addition  to  H2/NO2  mixtures.  The 
top  figure  shows  the  addition  of  10. 20.  30. 40.  50,  and  250  ppm  O2  to  the  initial  condi¬ 
tions  of  EXP  4  of  TABLE  6.  The  results  indicate  O2  addition  to  have  an  insigmficant  af¬ 
fect  on  the  rate  of  NO2  consumption.  The  initial  temperature  of  the  system  is  870  K  which 
at  3  atmospheres  is  above  the  extension  of  the  second  explosion  hmit  of  H2/O2  mixtures. 
As  a  result,  the  rate  of  H  +  O2  +  M  exceeds  that  of  H  -i-  O2  and  significant  branching  from 
the  addition  of  O2  is  prevented.  At  3  atmospheres,  the  transition  temperature  where 
branching  in  the  system  overtakes  stabilization  of  the  hydroperoxyl  radical  is  approxi¬ 
mately  970  K.  The  bottom  half  of  FIGURE  40  shows  the  kinetics  of  the  same  mixtures, 
except  for  an  initial  temperature  of  1070  K.  which  is  below  the  extension  of  the  second 
limit.  Note  that  in  this  case,  O2  does  not  have  a  significant  effect  due  to  the  rate  of  reaction 
NO2  +  H  exceeding  that  of  the  reaction  H  -t-  O2. 

Nitric  dioxide  decomposition  was  also  smdied  over  a  limited  set  of  conditions  (T  >  950 
K).  Due  to  complications  with  the  mixer-diffuser  design,  comparisons  of  these  data  with 
model  predictions  have  not  been  earned  out.  The  present  mixer/diffuser  design  incorpo¬ 
rates  some  stainless  steel  surfaces  which  apparently  cause  N02  surface  reactions  at  these 
conditions.  With  the  new  quartz  mixer-diffuser  section  in  place,  these  experiments  will  be 
repeated.  Modeling  studies  show  the  decomposition  rate  to  be  controlled  by  three  reac¬ 
tions. 

NO  2  +  (9  ^  NO  +  0  +  !W 
NO  2  ^M<r^NO  +  O^M 
N  O  2"^  N  O  2  ^  0  O . 

S3  Discussion 

Although  the  reactions  of  H2/NO2  system  are  generally  well  understand,  their  rates  in  the 
temperature  and  pressure  ranges  of  the  systems  studied  here  appear  to  need  further  refine¬ 
ment.  The  kinetics  of  the  present  systems  were  found  to  be  sensitive  to  reactions  which 
form  and  consume  HONO.  In  particular,  current  values  for  the  initiation  reacoon  may  be 
as  much  as  a  factor  of  5  to  large.  Because  this  system  forms  the  building  block  on  which 
other  mechanisms  for  the  decomposition  of  mtramines  are  to  be  constructed,  an  accurate 
desenpuon  of  the  kinetics  is  important. 
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FIGURE  39.  Calculated  temperature  and  mass  fraction  profiles  for 
the  reaction  of  H;  and  NOi  for  the  conditions  of  EXP  4  of  TABLE  6 
with  reaction  28  decreased  by  a  factor  of  5 


6.0  FUTURE  WORK 

Research  will  continue  on  the  H2/N2O  and  H2/NO2  systems.  Additional  fuels  which  re¬ 
main  to  be  studied  include  CO.  CH2O.  CH4.  C2H4.  HCN.  and  NH3.  We  also  expect  to 
have  in-situ  optical  diagnosuc  capability  for  OH  radicals  operational  in  the  variable  pres¬ 
sure  flow  reactor  during  the  next  year.  The  kinetic  models  for  H2/N2O/NO2  systems  will 
be  completed  and  the  expenmental  and  modeling  of  the  kinetics  on  the  additional  fuels 
listed  above  will  be  initiated. 


FIGURE  40.  Calculated  NOi  mass  fraction  profiles  for  the  initial 
conditions  of  EXP  4  of  TABLE  6  with  various  quantities  of  Oi 
added  at  870  K  (top)  and  1070  K  (bottom). 
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8.0  APPENDIX 


A.O  The  Atmospheric  Pressure  Flow  Reactor  Experiment 

An  atmospheric  pressure  flow  reactor  used  to  obtain  the  intermediate  temperature,  atmo¬ 
spheric  pressure  kinetic  data  is  also  available  for  reaction  studies.  This  reactor  is  designed 
to  approximate  a  steady,  one-dimensional  and  adiabatic  reacting  flow  with  negligible  radi¬ 
al  gradients  of  species  and  energy.  Gaseous  samples  are  withdrawn  from  the  reactor  at  up 
to  fifteen  discrete  locations  along  its  longitudinal  axis  with  a  variable  temperature,  convec¬ 
tion  quenched,  isothermal  sample  acquisition  and  storage  system.  The  sampling  system 
can  operate  at  temperatures  between  25'-^C  and  80°C.  Elevated  storage  temperature  pre¬ 
vents  the  condensation  of  high  boiling  point  intermediate  species  and  prevents  polymer¬ 
ization  reactions.  The  storage  system  interfaces  with  a  microprocessor  controlled  gas 
chromatograph  for  fully  automatic  analysis  of  the  stored  samples.  In  addition,  a  continu¬ 
ous  analysis  train  has  been  developed  for  determining  CO.  CO2.  O2  and  H2  concentrations 
in  the  sample  stream  from  which  the  discrete  samples  are  obtained.  The  flow  reactor,  its 
sampling  system,  and  the  analytical  instrumentauon  are  discussed  in  detail  below. 

A.l  The  Flow  Reactor 

The  atmosphenc  pressure  flow  reactor,  developed  in  this  laboratory  and  extensively  docu¬ 
mented  elsewhere  [1  -  3]  is  shown  schematically  in  FIGURE  1.  Fundamentally,  it  consists 
of  a  cylindrical  quartz  duct  through  which  hot  inen  carrier  gas  (nitrogen)  flows  at  high  ve¬ 
locities.  A  small  fraction  of  the  carrier  (S-lO'^c)  is  divened  through  an  arc  plasma  dis¬ 
charge.  which  heats  and  panially  lomzes  it.  The  hot,  panially  ionized  gas  is  premixed  with 
the  bulk  of  the  nitrogen  in  a  plenum  chamber,  where  recombination  results  in  metastable 
mtrogen  molecules  (.N'2*).  Collisional  de-excitation  of  the  metastable  .N2  is  necessary  to 
prevent  its  funher  reaction  with  O2: 


N2*  +  O2  —  2  NO 

and  to  promote  this  de-excitation  the  carrier  flow  continues  through  a  high  surface  area, 
ceramic  honeycomb  section,  and  emerges  as  men  N2.  Because  the  ceramic  section  has  a 
large  thermal  capacity,  the  temperature  of  the  nitrogen  leaving  the  section  is  constant.  (1 
K.  relative)  even  though  fluctuations  in  the  power  output  from  the  plasma  torch  itself  (2%) 
results  in  temperature  fluctuations  in  the  nitrogen  entering  the  section.  Desired  reaction 
temperatures  (800-1300  K)  are  achieved  by  vanation  of  the  split  ratio  of  the  nitrogen  flow 
and  of  the  power  to  the  plasma  torch. 

Oxygen  is  injected  into  the  earner  gas  upstream  of  a  converging/ diverging  inlet  section 
and  IS  well  mixed  with  the  carrier  gas  before  it  reaches  the  throat  of  the  section.  At  that 
point,  gaseous  fuel  is  injected  through  four  radially  opposed  quartz  injector  tubes.  The 
techruque  used  to  pre-vaporize  liquid  fuel  species  is  desenbed  later. 

The  high  velocity  radial  injection  of  gaseous  fuel  into  the  flowing  earner  provides  for  rap¬ 
id  mixing  of  the  fuel  with  the  oxygen/carrier  gas  stream.  Through  proper  control  of  carrier 


gas  temperature,  mixture  ratio,  and  total  gas  flow  rate,  a  steady  reacuon  zone  extending 
over  approximately  a  meter  in  length  can  be  established  which  represents  total  chemical 
reaction  times  between  10  and  150  ms.  Mean  velocity  along  the  axis  of  the  reactor  duct  is 
expenmentally  determined  using  static  pilot  probes  [4]  in  reacting  flows  and  hot  wire  ane- 
mometry  in  non-reacting  flows.  The  use  of  small  imtial  concentrations  of  both  fuel  and  ox¬ 
idizer  and  the  spreading  of  the  reaction  zone  over  a  meter  gives  rise  to  longitudinal 
gradients  small  enough  for  diffusion  of  mass  and  energy  along  the  length  of  the  duct  to  be 
neglected  with  respect  to  convection  [5].  Furthermore,  the  dilute  nature  of  the  reaction  and 
in  particular,  the  limitations  therefore  imposed  on  the  turbulent  fluctuations  of  temperature 
from  the  mean  (1%)  limits  any  turbulent  fluctuation  coupling  [6]  with  chemical  kinetics  to 
be  negligible.  Measured  kinetics  have  been  found  to  be  independent  of  flow  Reynolds 
number  within  the  limits  of  experimental  accuracy  (10%).  The  relatively  large  duct  diam¬ 
eter  ( 10. 16  cm)  results  in  a  small  surface  to  volume  ratio  and  minimizes  wall  effects. 

Adiabaticity  is  approximated  by  heating  and  insulating  the  inlet  section  and  main  reactor 
duct,  so  that  the  walls  are  quickly  equilibrated  to  the  inlet  temperature  of  the  earner  gas. 
The  waU  temperature  is  monitored  at  four  locations  along  the  length  of  the  duct  by  Pt/Pt 
13%  Rh  thermocouples  placed  in  an  air  gap  between  the  heaters  and  the  duct. 

The  flow  reactor  experiment  requires  that  fuels  be  in  a  gaseous  state,  and  thus  any  liquid 
material  is  pre-vaponzed  using  an  evaporator.  The  evaporator  is  shown  schematically  in 
FIOLTE  2.  Nitrogen  pressure  on  the  fuel  in  the  storage  tank  forces  metered  amounts  of  liq¬ 
uid  fuel  through  an  orifice  onto  a  heated  porous  cylinder,  which  is  at  a  temperature  above 
the  boiling  point  of  the  fuel.  The  porous  cylinder  vaporizes  the  fuel  which  is  simultaneous¬ 
ly  mixed  with  a  small  amount  of  heated  nitrogen  in  the  evaporator  and  subsequently  di¬ 
rected  to  a  manifold  which,  in  turn,  feeds  the  radial  injector  rubes.  Eventually,  a  steady 
state  is  established  in  which  the  liquid  fuel  flow  into  the  evaporator  is  equal  to  the  gaseous 
fuel  flow  out.  In  order  to  prevent  condensation  in  the  transfer  lines  leading  into  the  reactor, 
the  flow  lines  from  the  evaporator  are  heated  and  insulated. 

In  summary,  the  flow  reactor  provides  a  well-defined  and  well-monitored  reaction  zone 
which  can  be  accurately  represented  as: 

•  constant  pressure  (small  energy  release,  low  Mach  number  flow  system) 

•  one-dimensional 

•  quasi-steady 

•  homogeneous  (wall  effects  negligible) 

•  adiabatic 

•  with  negligible  longitudinal  diffusion  of  mass  and  energy 

These  operation  charactenstics  and  the  reproducibility  of  the  experimental  data  will  be 
discussed  later  in  this  appendix, 

.4.2  Data  .Acquisition 

.\.2.1  Extractive  and  Immersion  Measurements 


Mean  temperature,  and  species  profiles  are  measured  along  the  length  of  the  reactor  duct. 
The  temperature  is  measured,  with  a  silica  coated  Pt  30%  Rh/Pt  6%  Rh  alloy  thermocou¬ 
ple  with  a  bead  diameter  of  38  micrometers  and  a  spatial  resolution  of  about  twenty  bead 
diameters.  Dryer  [2]  has  discussed  this  technique  at  length  in  his  Ph.D.  thesis. 

The  gas  sampling  techniques,  in  pnnciple,  are  also  those  used  in  earlier  flow  reactor  stud¬ 
ies:  combustion  gases  are  withdtawn  continuously  from  the  reactor  with  a  hot  water 
cooled,  water  convection  quenched  probe.  The  sample  stream  can  be  analyzed  continu¬ 
ously  and  discrete  volumes  can  also  be  extracted  from  the  sample  stream  and  stored  for 
later  analysis  by  a  gas  chromatograph. 

The  gas  sample  probe  itself  appears  in  Section  A  of  RGURE  1  and  is  constructed  of  three 
concentric  type  316  stainless  steel  tubes.  The  two  outer  tubes  provide  passages  for  the  hot 
water  (60-90'<r)  used  to  convection  cool  the  flowing  gas  sample  and  maintain  the  mechan¬ 
ical  integrating  of  the  probe.  The  gas  sample  quench  is  accomplished  purely  by  convection 
heat  transfer  at  a  rate  greater  than  10'  Kh  (estimated:  see  [2]).  The  sample  is  transferred 
from  the  probe  to  the  sampling  assembly  through  heated  (thermostated)  teflon  lines  by  a 
heated  stainless  steel  bellows  pump.  The  flow  exiting  the  pump  is  at  approximately  2  at¬ 
mospheres  pressure  supplying  sample  flow  to  one  or  several  multiloop  sampling  valves 
(described  later)  and  the  continuous  analysis  tram. 

The  continuous  analysis  system  (FIGURE  3)  provides  two  principle  advantages  in  conduct¬ 
ing  flow  reactor  research: 

•  It  permits  accurate  placement  of  the  reaction  zone  within  the  reactor  and  thus  optimize 
the  experiment  for  the  ranges  of  extent  of  reaction  to  be  studied. 

•  It  relieves  several  complications  imposed  when  simultaneous  gas  chromatographic 
analysis  of  light  permanent  gases  H2.  CO.  CO2  and  higher  molecular  weight  hydrocar¬ 
bons  IS  required. 

Carbon  monoxide  and  carbon  dioxide  arc  determined  using  non-dispersive  infrared  analy¬ 
sis  (Beckman.  Inc.).  Oxygen  is  determined  using  a  thermomagnetic  techmque  (Leeds  and 
Northrup).  A  measurement  of  hydrogen  is  obtained  by  catalytic  conversion  of  the  dned 
sample  stream  and  determination  of  the  resulting  water  vapor  by  a  Dupont  Acoustic  water 
vapor  analyzer. 

What  has  recently  advanced  the  analytical  technology  to  deal  with  the  oxidation  of  higher 
carbon  number  paraffins  and  aromatics  is  the  development  of  a  system  which  would  si¬ 
multaneously  acquire  and  store  up  to  two  discrete,  heated  samples  at  each  longitudinal  po¬ 
sition  sampled  m  the  flow  reactor.  This  technique  permits  simultaneous  use  of  a  Hewlett 
Packard  HP5840  gas  chromatograph  (GC)  for  quantitative  analysis  of  the  sample  for  hy¬ 
drocarbons  and  of  a  HP5895  gas  chromatography/mass  spectrometer  (GCAIS)  for  the  re¬ 
spective  specie  identifications  (both  by  column  retention  time  and  mass  fractionation 
patterns). 

Successful  discrete  chemical  sampling  involves  satisfactory  accomplishment  of  each  of 
the  following  goals: 
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•  isolation  of  a  flow  of  reactive  media  from  a  specific  point  in  the  reaction  environment. 

•  prevention  of  further  chemical  or  thermodynamic  alteration  of  the  contained  species. 

•  introduction  of  a  defined  quantity  of  this  sample,  intact,  to  an  analytical  instrument. 

The  heated  sampling  system  [3]  which  is  now  commonly  used  has  several  distinct  advan¬ 
tages.  It  interfaces  with  all  analytical  equipment  in  a  manner  completely  analogous  to  its 
GC  interface,  and  allows  ready  identification  of  species  which  are  impractical  to  identify 
with  a  GC  alone  and  which  would  be  lost  to  condensation  of  room  temperature. 

A  heated  storage  temperature  of  70°C  has  a  second  advantage  related  to  the  integrity  of  the 
stored  sample:  by  preventing  polymerizafcn  reactions  associated  with  the  storage  of 
formaldehyde  at  room  temperature,  it  should  allow  accurate  formaldehyde  analysis.  Con¬ 
densation  of  other,  low  vapor  pressure  species  is  also  prevented  by  heated  storage. 

Funhermore.  the  heated  sampling  system  has  the  advantages  of  preventing  exposure  of 
stored  samples  to  light  and  to  adsorptive  surfaces  (e.g.  mercury,  synthetic  rubber),  both  of 
which  occurred  in  the  older  system.  All  surfaces  which  contact  sample  gases  in  the  heated 
sampling  system  are  stainless  steel,  so  the  exposure  to  light  that  resulted  from  storage  in 
glass  bottles  on  an  exposed  sampling  rack  is  circumvented.  Furthermore,  gases  are  pres¬ 
surized  for  injection  with  a  stainless  steel  welded  bellows  pump  and  not  with  the  mercury 
system  previously  in  use. 

Operation  of  the  heated  sampling  storage  and  chromatographic  injection  facilities  are  as 
foUows  (FIGURE  4).  Combustion  gases  are  withdrawn  from  the  reactor  through  the  probe 
and  a  length  of  heated  line  by  a  heated  stainless  steel  welded  bellows  pump.  The  bellows 
pump  raises  the  pressure  of  the  sampled  gases  to  approximately  two  atmospheres  pressure, 
and  the  flow  proceeds  through  a  heated  7  micro  filter  and  temperature  equilibrium  coils  to 
two  sample  storage  vales.  The  temperature  equilibration  coil  and  the  storage  valves  are 
contained  in  separate,  portable,  thermostated  oven  assemblies. 

Each  storage  valve  itself  is  a  Valeo  32  port.  16  position  pneumatically  actuated  valve  fined 
with  sixteen  stainless  steel  sample  storage  loops.  Each  of  the  valve  loops  has  a  pre-deter- 
mined  volume  of  12  cc.  When  the  valve  is  actuated,  the  flow  is  switched  from  one  loop  to 
the  next,  a  defined  quantity  of  sample  is  trapped  in  the  sealed  loop  (P.  V  and  T  are  all  de¬ 
fined),  and  the  next  sample  loop  is  purged.  Up  to  fifteen  discrete  samples  can  be  stored  in 
this  manner,  with  the  remaimng  loop  open  to  the  flow  path.  Once  filled,  the  multiposition 
valve  system  is  carried  (m  its  oven)  to  a  gas  chromatograph,  where  the  loops  are  flushed, 
one  by  one,  into  a  chromatographic  gas  sample  injection  system  (RGURE  5)  and  injected 
onto  the  GC  columns  for  analysis. 

A  heated  1  cc  single  loop  sample  injection  valve  is  used  to  determine  the  sample  size  for 
actual  analysis.  This  sample  injection  valve  is  connected  on  the  inlet  side  to  the  multiposi- 
tion  valve  (i.e.  to  the  stored  samples)  through  a  length  of  heated  line,  and  on  the  outlet  side 
to  a  manostatAacuum  system.  When  the  multiposition  valve  is  actuated  and  a  new  storage 
loop  is  exposed  to  the  sample  injection  system,  the  gases  in  the  12  cc  loop,  at  2  atmo¬ 
spheres,  flow  through  the  1  cc  secondary  sample  injection  loop  until  the  pressure  in  the 
system  is  equilibrated  to  0.5  atmospheres  (controlled  by  a  cartesian  manostai).  Sufficient 


excess  sample  is  stored  in  the  multiposition  valve  loops  to  assure  complete  displacement 
of  any  residual  gases  in  the  secondary  sample  valve  and  transfer  lines.  The  0.5  atm.  sam¬ 
ple  in  the  secondary  sample  valve  is  then  injected  onto  the  gas  chromatograph  columns  for 
analysis.  Actuation  of  both  the  multiposition  valve  and  the  secondary  sample  valve  is  con¬ 
trolled  by  a  microprocessor  indigenous  to  the  gas  chromatograph,  so  fully  automatic  anal¬ 
ysis  of  successive  samples  is  possible.  This  procedure  not  only  assures  a  constant 
systematic  error,  but  minimizes  sample  storage  times  and  reduces  dramatically  the  man¬ 
power  which  must  be  devoted  to  the  analysis  of  samples.  Previous  (manual)  injection 
techniques  required  operator  attention  throughout  the  25+  hour  analysis  time  required  for 
a  complete  flow  reactor  analysis.  With  the  present  system  interfaced  to  the  gas  chromato¬ 
graph  and  the  GC  microprocessor  preprogrammed,  no  additional  operator  attention  is  re¬ 
quired  for  complete  analysis. 

Because  the  sam.ple  injected  onto  the  GC  columns  is  set  in  the  same  loop  each  time,  the  ef¬ 
fects  of  loop  volume  vanability  are  eliminated.  Sample  storage  valve  leakage  problems 
are  also  circumvented;  since  samples  are  stored  at  a  positive  pressure,  gases  leak  out.  not 
in.  Any  loss  of  sample  at  the  small  rates  resulting  from  loop  to  atmosphere  and  loop  to 
loop  diffusion  is  not  reflected  in  the  sample  delivered  to  the  sampling  injection  valve, 
since  a  large  volume  of  sample  gas  is  initially  stored.  Cross  leak  distortion  of  the  stored 
sample  has  been  shown  to  be  less  than  1%.  even  for  the  lightest  gas  expected  in  hydrocar¬ 
bon  oxidation  studies  (CH).  Note  that  since  diffusion  rates  are  inversely  proportional  to 
the  square  root  of  the  molecular  weight  of  the  diffusing  species,  errors  for  most  species 
will  be  significantly  smaller. 

In  summary,  the  advantages  of  this  unique  sampling  storage  and  chromatographic  injec¬ 
tion  system  are; 

•  control  of  the  sample  acquisition  and  storage  temperatures  (up  to  about  80  C). 

•  sample  isolation  from  adsorptive  surfaces. 

•  sample  isolation  from  light  surfaces. 

•  automation. 

•  Additional  details  appear  in  the  recent  thesis  of  Euchner  [3]. 

A. 2.2  Reproducibility  of  Experimental  Data 

We  have  recently  performed  two  research  studies  [4.  7]  which  adequately  demonstrate  the 
nature  and  quahty  of  the  experimental  data  obtained  from  the  flow  reactor. 

In  a  study  on  the  elementary  combustion  chemistry  of  the  CO/H2/O2  system  [4].  we  have 
performed  extensive  experiments  which  involve  the  turbulent  flow  reactor  and  the  contin¬ 
uous  sampling  system.  In  the  performance  of  some  of  these  experiments  we  investigated 
the  ability  to  reproduce  the  same  experimental  initial  conditions  and  chemical  analysis 
data  by  attempting  to  repeat  the  same  expenment  on  size  occasions  over  a  penod  of  6 
months.  The  experimental  conditions  and  results  are  shown  in  FlGLTlE  6.  The  precision  of 
these  data  is  clearly  evident  and  indicates  both  the  ability  to  reproduce  experimental  van- 
ables  and  also  the  continuous  sampling  instrumentation  calibration  effects 


FIGURE  7  presents  the  comparison  of  elementary  kjnetics  calculations  performed  using  ac¬ 
cepted  literature  rate  constant  and  mechanism  with  experimental  measurements  on  the 
CO/H2/O2  system  [4]  (see  TABLES  5  and  6).  In  the  four  flow  reactor  experiments  (point 
data),  only  the  initial  reaction  temperature  was  changed.  The  three  calculations  (line  data) 
performed  show  very  good  agreement  with  the  experimental  data,  even  for  reaction  times 
approaching  0.6  seconds!  It  should  be  emphasized  that  there  was  no  performance  of  “em¬ 
pirical  fitting"  of  the  model  to  the  data,  and  no  effects  of  diffusive  and  turbulent  transport 
or  coupling  with  chemistry  were  considered  in  the  model.  These  comparisons  therefore 
support  the  “pure  chemistry”  interpretation  of  flow  reactor  data.  Finally  it  should  be  em¬ 
phasized  that  the  experimental  and  calculated  results  are  for  condition  crossing  the  explo¬ 
sion  limit  regime  for  CO/H2/O2  chemistry  at  one  atmosphere. 

A  second  set  of  data  which  demonstrate  the  precision  of  the  discrete  chemical  sampling 
and  gas  chromatographic  analytical  results  are  presented  in  FIGURE  8.  These  data  were 
performed  as  pan  of  a  recent  thesis  [7]  and  were  obtained  from  the  analysis  of  samples  (16 
samples  by  each  valve,  one  sample  in  each  valve  at  each  sampling  position)  taken  by  two 
sample  storage  valves  in  parallel  in  the  flow  reactor  sampling  systems.  The  16  loops  of 
each  sample  valve  were  then  analyzed  sequentially,  with  the  analysis  of  the  second  valve 
samples  following  the  first.  The  chemical  analysis  thus  spanned  a  period  of  some  72  hours 
after  the  flow  reactor  expenment  had  been  concluded.  The  reproducibility  of  this  data  is 
again  excellent. 

A.2.3  Optical  Measurements 

The  flow  reactor  has  optical  access  at  three  axial  positions  for  line  of  site  absorption  mea¬ 
surements.  Linteris  [8].  has  shown  that  OH  can  be  quantitatively  measured  in  the  CO/ 
H2O/O2  reactor  system  at  flow  reactor  temperatures  using  resonance  absorption,  and  that 
the  results  compare  favorably  with  those  predicted  [4].  In  order  to  obtain  measurements  at 
other  locations  in  the  reactor,  a  unique  optical  probe  was  developed  [8].  based  on  180'  la¬ 
ser  induced  fluorescence  (LIF).  The  techmque  has  been  shown  to  be  feasible  and  to  pro¬ 
duce  OH  profiles  which  could  be  successfully  calibrated  by  simultaneous  resonance 
absorption  measurements.  Examples  of  the  results  are  reported  in  RGUTIE  9. 

A.3  REFERE.NCES 

41. Crocco,  L.,  Classman.  I.  and  Smith.  I.E..  Jet  Propulsion.  1266  (1957). 

42. Dryer.  F.L..  Ph.D.  Thesis.  Department  of  Aerospace  and  Mechanical  Sciences.  Prince¬ 
ton  Umversity.  Princeton.  NJ  (1972). 

43. Euchner,  J.A..  M.S.E.  Thesis,  Mechanical  and  Aerospace  Engineering  Department. 
Princeton  University.  Princeton,  NJ  (1980). 

44. Yetter.  R.A.,  Ph.D.  Thesis,  Mechanical  and  Aerospace  Engineering  Department.  Princ¬ 
eton  University.  Princeton,  NJ  (1985). 

45.  Westbrook,  C.K..  Creighton.  J..  Lund,  C.  and  Dryer.  F.L.,  J.  Phys.  Chem..  81. 2542 
(1977). 

46. Gouldin.  F.C..  Comb.  Sci.  Tech..  9.  17  (1974). 


47. Proscia.  W.M..  M.S.E.  Thesis.  Mechanical  and  Aerospace  Engineering  Depanment. 
Princeton  University,  Princeton,  NJ  (1983). 

48. Linteris,  G.T.,  Ph.D.  Thesis,  Depanment  of  Mechanical  and  Aerospace  Engineering 
Department,  Princeton  University,  Princeton,  NJ  (1990). 


A.4  TABLES 

TABLE  7.  €0112/02  Thermochemical  Parameters* 


SPECIES 

4H/.298 

Sms 

Opjoo 

Cp.400 

^p.500 

Cp.600 

Cp.800 

Cp.lOOO  f^p.isoo 

H 

52.10 

27.39 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

4.97 

0 

59.56 

38.47 

5.23 

5.14 

5.08 

5.05 

5.02 

5.00 

4.98 

OH 

9J2 

43.88 

7.15 

7.10 

7.07 

7.06 

7.13 

7.33 

7.8 

H2 

0.00 

31.21 

6.90 

6.96 

7.00 

7.02 

7.07 

7.21 

7.73 

O’ 

0.00 

49.01 

7.01 

7.22 

7.44 

7.65 

8.07 

8.35 

8.72 

H2O 

-57.80 

45.10 

8.00 

8.23 

8.44 

8.67 

9.22 

9.87 

U.26 

HO, 

3.00 

54.43 

8.36 

8.95 

9.48 

9.96 

10.75 

11.37 

1234 

HjO, 

-32.53 

55.66 

10.41 

11.44 

12.34 

13.U 

14.29 

15.21 

16.85 

CO 

-26.42 

47.21 

6.95 

7.03 

7.14 

7.27 

7.61 

7.95 

8.41 

0 

0 

1,^ 

-94.06 

51.08 

8.91 

9.86 

10.65 

11.31 

12.32 

12.99 

13.93 

HCO 

10.40 

53.66 

8.24 

8.78 

9.28 

9.77 

10.74 

11.52 

1236 

N, 

0.00 

45.77 

6.95 

7.01 

7.08 

7.19 

7.50 

7.83 

8.32 

“  units  are  kcaVmol  for  AHy  29s  callmol-K  for  S  and  Cp 


TABLE  8.  CO  H2'02  Reaction  Mechanism 


^^^298 


H2-O2  Chain  Reactions 

1,  H+O2  =  O+OH  16.77 

2.0+H2=H+0H  1.86 

3.  OH+H2  =  H+H2O  -15.02 

4.  OH+OH  =  O+H2O  16.87 

H2-O2  Dissociation/Tiecombination  Reactions 

5.  H2+M  =  H-rH+.M  104.1 

6.  O+O+M  =  O2+M  -119.1 

7.  O-^H+M  =  OH-.^M  -102J 

8.  H+OH+M  =  H.O-kM  -U9.2 

Formation  and  Consumption  of  HO2 

9.  H-kO,+M  =  HO,+M  -49.07 


log(Af) 

"r 

Ea.r 

L'F 

fRANGE 

14.28 

0.00 

16.44 

2 

962-2577 

4.71 

2.67 

6.29 

1.5 

297-2495 

8.33 

1.51 

3.43 

1.5 

250-2581 

k=5.46xl0“ 

exp(0.00149xT)' 

2.5 

250-2000 

19.66 

-1.40 

104.38 

3.0 

600-2000 

15.79 

-0.50 

0.00 

1.3 

2000-10000 

18.67 

-1.00 

0.00 

10.0 

22.35 

-2.00 

0.00 

2. 

1000-3000 

19.83 

-1.42 

0.00 

3. 

200-2000 

10.  HO2+H  =  H.+O, 

-55.07 

13.82 

0.00 

2.13 

2. 

298-773 

U.  HO2+H  =  OH+OH 

-36.44 

14.23 

0.00 

0.87 

2. 

298-773 

12.  H0,+0  =  OH+O, 

-53.22 

13.24 

0.00 

-0.40 

1.2 

200-400 

13.  HO2+OH  =  H2O+O2 

-70.09 

16.16 

-1.00 

0.00 

2. 

298-1400 

Formation  and  Consumption  of  H2O2 

14.  HO2+HO2  =  H2O2+O2 

-38.51 

12.48 

0.00 

1.39 

3. 

650-800 

15.  H2O2+M  =  OH+OH+M 

51.14 

17.08 

0.00 

45.50 

2. 

700-1500 

16.  H2O2+H  =  HjO+OH 

-68.02 

13.00 

0.00 

3.59 

3. 

283-800 

17.  H2O2+H  =  H2+HO2 

-17.56 

13.68 

0.00 

7.95 

5. 

250-800 

19.  HjOj+OH  =  H2O+HO2 

-32.57 

12.85 

0.00 

1.43 

2. 

98-800 

Oxidation  of  CO 

20.  CO+O+M  =  CO2+M 

127.01 

13.40 

0.00 

-4.54 

4. 

1000-3000 

21.  CO-rOj  =  CO2  +  O 

-8.08 

12.40 

0.00 

47.69 

2. 

500-3000 

22.  CO+OH  =  CO2+H 

-24.85 

k=6.75xl0’®exp[0.000907xT]‘' 

1.5 

250-2500 

23.  CO+HO2  =  CO2+OH 

-60.79 

13.78 

0.00 

22.95 

3. 

700-1000 

Formation  and  Consumption  of  HCO 

24.  HCO+M  =  H+CO+M 

15.27 

17.27 

-1.0 

17.00 

5. 

637-832 

25.  HCO+O,  =  CO+HO, 

-34.30 

12.62 

0.00 

0.00 

5. 

295-713 

26.  HCO+H  =  CO+H, 

-88.87 

13.86 

0.00 

0.00 

2. 

296-418 

27.  HCO+O  =  CO+OH 

-87.02 

13.48 

0.00 

0.00 

3. 

28.  HCO+OH  =  CO+HjO 

103.09 

13.48 

0.00 

0.00 

3. 

... 

Sinits  are  cm^-moU-s-kcal,  k  =  AT"  exp(-Eg/RT) 

^  [M]  =  [N2]+[Hl+[0]+[0H]+2.5[H2l+[02]+12[H:0]+lH02]+[Hp2]+1.9[C0]+3.8[C02]+[HC0] 
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A.5  FIGURES 


Inlet  Section  Test  Section 


FIGURE  1.  Schematic  diagram  of  flow  reactor  and  flow  measurements,  respectively.  T,  P  and  F 
indicate  temperature,  pressure 
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FIGURE  2.  Evaporator  schematic 


FIGURE  3.  The  continuous  analysis  system 


sample  acquisition  system  schematic 


FIGURE  4.  Sample  acquisition  system  schematic 
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FIGURE  7.  Mean  CO  concentration  profiles  for  COHiO  Oj  mixtures  with  the  initial  temperature 
ranging  from  917  to  952K.  Initial  conditions;  circles;  X(CO)  =  0.00895,  X(Ot)  =  0.00438,  X(HiO)  = 
0.01235,  X(Nt)  =  0.97432,  T  =  952K,  p  =  1  atm,  m  =  6.39  gm  s.  squares;  X(CO)  =  0.00905,  X{(5^  = 
0.00442.  XiHiHjO)  =  0.01324,  XiN^)  =  0.97329,  T  =  943K.  p  =  1  atm,  m  =  6.26  gm/s,  diamonds;  X(CO) 
=  0.009i5,  XfO^  =  0.00448,  X(H,0)  =  0.0133.  X(N,)  =  0.97307,  T  =  932K:,  p  =  1  atm,  m  =  6.28  gm's, 
triangles;  X(CO)  =  0.0093,  X(02)’  =  0.00455,  XIH^O)  =  0.0147,  XIN,)  =  0.97145,  T  =  917K,  p  =  1  atm, 
m  =  6.17  gm/s. 
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